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ABSTRACT 

We attempt to get a physical insight into grain ahgnment processes by studying basic 
properties of radiative torques (RATs). For this purpose we consider a simple toy 
model of a helical grain that reproduces well the basic features of RATs. The model 
grain consists of a spheroidal body with a mirror attached at an angle to it. Being 
very simple, the model allows analytical description of RATs that act upon it. We 
show a good correspondence of RATs obtained for this model and those of irregular 
grains calculated by DDSCAT. Our analysis of the role of different torque components 
for grain alignment reveals that one of the three RAT components does not affect the 
alignment, but induces only for grain precession. The other two components provide 
a generic alignment with grain long axes perpendicular to the radiation direction, if 
the radiation dominates the grain precession, and perpendicular to magnetic field, 
otherwise. The latter coincides with the famous predictions of the Davis-Greenstein 
process, but our model does not invoke paramagnetic relaxation. In fact, we identify a 
narrow range of angles between the radiation beam and the magnetic field, for which 
the alignment is opposite to the Davis-Greenstein predictions. This range is likely to 
vanish, however, in the presence of thermal wobbling of grains. In addition, we find 
that a substantial part of grains subjected to RATs gets aligned with low angular 
momentum, which testifies, that most of the grains in diffuse interstellar medium do 
not rotate fast, i.e. rotate with thermal or even sub-thermal velocities. This tendency 
of RATs to decrease grain angular velocity as a result of the RAT alignment decreases 
the degree of polarization, by decreasing the degree of internal alignment, i.e. the 
alignment of angular momentum with the grain axes. For the radiation-dominated 
environments, we find that the alignment can take place on the time scale much 
shorter than the time of gaseous damping of grain rotation. This effect makes grains 
a more reliable tracer of magnetic fields. In addition, we study a self-similar scaling of 
RATs as a function of X/aeff. We show that the self-similarity is useful for studying 
grain alignment by a broad spectrum of radiation, i.e. interstellar radiation field. 

Key words: ISM- Magnetic fields- polarization, ISM: dust-extinction 



1 INTRODUCTION 

Magnetic fields play a crucial role in many astrophysical 
processes, e.g. star formation, accretion of matter, trans- 
port processes, including heat conduction and propagation 
of cosmic rays. One of the easiest ways to study magnetic 
field topology is via polarization of radiation arising from 
extinction or/and emission by aligned dust grains. The new 
instruments like Scuba II (Bastien, Jenness & Molnar 2005), 
Share II (Novak et al. 2004) and an intended polarimeter for 
SOFIA open new horizons for tracing of astrophysical mag- 
netic fields with aligned grains. 

In this situation it is unacceptable that the processes 
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of grain alignment are not completely understood (see re- 
view by Lazarian 2003). The enigma that surrounds grain 
alignment since its discovery in 1949 (Hall 1949; Hilner 1949) 
makes the interpretation of the polarization in terms of mag- 
netic fields somewhat unreliable. The failure of grains to 
align at high optical depths was discussed, for instance, in 
Goodman (1995). 

A recent progress in understanding of the grain align- 
ment physics removed many questions, but have not reme- 
died the situation completely. Among the milestones let us 
mention the recent revival of interest to radiative torques 
(henceforth RATs) . Introduced by Dolginov & Mytrophanov 
(1976), those torques, that arise from the interaction of ir- 
regular grains with a flow of photons, were essentially forgot- 
ten till Draine & Weingarter (1996, 1997, henceforth DW96, 
DW97, respectively) provided quantitative numerical stud- 
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ies. While Dolginov & Mytrophanov (1976) were somewhat 
vague on what makes RATs important for a grain, DW96 
demonstrated that their arbitrary chosen irregular grains 
exhibit dynamically important RATs when subjected to a 
typical interstellar radiation field (ISRF) . Very importantly, 
Bruce Draine incorporated RATs into the publicly available 
DDSCAT code (Draine & Flatau 1994), which stimulated 
a further progress in the field. First laboratory studies of 
RATs were reported in Abbas et al. (2004). 

The renewed interest to RATs coincided with a crisis 
of the paramagnetic alignment as it is described in text- 
books (Purcell 1979; Spitzer & McGlynn 1979; Mathis 1986). 
Lazarian & Draine (1999a) (hereafter LD99a) identified new 
elements of grain dynamics, which they termed "thermal 
hipping" and "thermal trapping" . Due to thermal wobbling 
arising from the dissipative coupling of grain vibrational and 
rotational degrees of freedom (Lazarian 1994; Lazarian & 
Roberge 1997) grains smaller than a critical radius Oc flip fre- 
quently and thus average out uncompensated torques. These 
torques, that were first discussed by Purcell (1979), were 
considered essential to make otherwise inefficient paramag- 
netic alignment (Davis & Greenstein 1951; Jones & Spitzer 
1967) to account for the polarimetric observations. A new 
dissipative coupling mechanism related to nuclear spins of 
constituent atoms, that was described in Lazarian & Draine 
(1999b) (henceforth LD99b), resulted in ac larger than the 
typical cut-off scale for grains in diffuse interstellar medium 
(ISM). As the other mechanisms, e.g. mechanical alignment 
(Gold 1951, review by Lazarian 2003 and references therein), 
have their limitations, this made the RAT alignment the 
only viable mechanism to explain the ubiquity of interstellar 
polarization (DW97) and possibly polarization arising from 
aligned dust in other astrophysical environments (Lazarian 
2007). 

The successes of RATs include a more recent work by 
Cho & Lazarian (2005), where a substantial increase of the 
RAT efficiency with the grain size was established. This work 
explained the sub-millimeter polarization data for quiescent 
starless cores (Wart-Thompson et al. 2000) by appealing to 
the differential RAT alignment of large grains. For such cores 
the analysis of all other mechanisms in Lazarian, Goodman 
& Myers (1997) predicted only marginal degrees of align- 
ment. The studies elaborating the approach in Cho & Lazar- 
ian (2005), e.g. Pelkonen, Juvela & Padoan (2007), BetheU et 
al. (2007) , provided theory-motivated predictions of the de- 
gree of alignment for numerically simulated molecular clouds 
and cores. 

However, the above explanation as well as other expla- 
nations (see Lazarian 2003) are based on the plausibility of 
arguments, rather than on the rigorous RAT alignment the- 
ory. Indeed, DW96 considered RATs as a means of spin-up. 
This induced a naive explanation of RAT alignment action 
that could be perceived in some of the papers that followed 
the DW96 study. There it was assumed that RATs were 
proxies of the Purcell's torques (1979), that arise from the 
action of photons, rather than from the action of II2 forma- 
tion over catalytic sites, as in the original mechanism. While 
the Purcell's torques depend on the resurfacing and therefore 
short lived, RATs depend on grain shape and can be long- 
lived. As a result, long-lived fast rotation of paramagnetic 
grains should induce good paramagnetic alignment (Purcell 



1979). This understanding of RATs is not correct, as it is 
clear from a more careful reading of DW96 and DW97. 

In fact, RATs can be subdivided into the parts that 
arise from isotropic and anisotropic radiation fluxes. The 
part arising from anisotropic radiation, for which we adopt 
a shorthand notation "isotropic part", is, indeed, similar to 
the Purcell's torques. The "anisotropic part" is, however, 
both usually stronger and has properties different from the 
Purcell's torques. The major difference arises from the fact 
that RATs are defined in the laboratory, rather than in the 
grain coordinate system. Thus, the presence of even a small 
anisotropic component of radiation, which is a natural con- 
dition for any realistic astrophysical system, is bound to 
change the dynamics of grain. Note, that the alignment by 
the anisotropic radiation was first discussed by Doginov & 
Mitraphanov (1976). They, however, concluded that prolate 
and oblate grains can be aligned differently. Lazarian (1995) 
took into account internal relaxation and claimed that both 
prolate and oblate grains should be aligned with longer axes 
perpendicular to magnetic field. Nevertheless, the theory 
lacked a proper description of RATs. 

DW97 demonstrated numerically that in the presence 
of anisotropic radiation the grains can be aligned by RATs 
in respect to magnetic field on the time scales much shorter 
than the time scale for paramagnetic alignment In gen- 
eral, the magnetic field for the RAT alignment acts through 
inducing fast Larmor precession; the alignment potentially 
may happen both with long grain axes parallel and perpen- 
dicular to magnetic field. Only the latter is consistent with 
polarimetric observations, however (see Serkowski, Mathew- 
son & Ford 1975). 

In the DW97 study, the alignment with longer grain 
axes perpendicular to magnetic field ("right alignment") 
happened more frequently than the grain alignment with 
longer grain axes parallel to magnetic field ("wrong align- 
ment"). This experimental evidence, based on a limited sam- 
pling, raised worrisome questions. Is this a general property 
of radiative torques or just a coincidence? Do we expect 
to see more of "wrong alignment" if the grain environment 
is different from the interstellar one? What are the chances 
that we are fooled by the "wrong alignment" while interpret- 
ing the polarimetry measurements in terms of the underlying 
magnetic fields? It seems necessary to address these ques- 
tions if interpreting polarimetry data in terms of underlying 
magnetic fields is sought. 

Analytical calculations played an important role for for- 
mulating the models of both paramagnetic and mechanical 
alignment (see Davis & Greenstein 1951, Jones & Spitzer 
1967, PurceU 1979, Spitzer & McGlynn 1979). Ahhough 
such calculations dealt with intentionally idealized models of 
grains, they allowed deep insight into the relevant physics. 
Dolginov & Mytrophanov (1976) attempted an analytical 
modeling for RATs. They used a model grain containing 
two ellipsoids connected together at an angle. The radiative 
torques were calculated for such a model grain by assuming 
that the wavelength is much larger than the grain size, i.e., 
in the Rayleigh-Hans approximation (Dolginov & Silantev 



^ DW97 identifies this time-scale with the gaseous damping time. 
In §5.5 we show that the alignment could happen much faster in 
the presence of strong radiation sources 
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1976). However, adopting their shape, we could not repro- 
duce numerically their analytical predictions for RATs. This 
induce us to seek analytical models that would correspond 
to the DDSCAT calculations. 

Our approach in the present paper is to provide a phys- 
ical insight into basic RAT properties. In §2 we explain why 
we consider only RATs due to anisotropic radiation. Then, 
we describe a simple grain model that is subject to RATs and 
allows analytical descriptions (^. In ^we present calcula- 
tions of RATs for a number irregular shapes and study the 
correspondence of their RATs with AMO. We also briefly 
consider possible generalizations of our model (^SJ. In !j6] 
we analyze the alignment for both this model and irregu- 
lar grains with respect to the radiation direction. In ^we 
study the alignment of our model and irregular grains with 
respect to magnetic field. Crossovers are studied in ^ while 
we identify the conditions for the magnetic field or the radi- 
ation direction to act as the axis of alignment in !j9l As, even 
with modern computers, the calculations of RATs for a vari- 
ety of wavelengths is time consuming, we address the ques- 
tion of the accuracy of presenting radiative torques as the 
function of the ratio of grain size to the wavelength X/a^ff 
in i jlOl The discussion of our results and the summary are 
provided in ijlll and §12, respectively. 



2 ISOTROPIC AND ANISOTROPIC RATS 

RATs can emerge even when the radiation field is isotropic. 
Devices similar to those used by Lebedev (1901) to measure 
radiation pressure experience torques in the presence of the 
isotropic radiation (see also the cartoon of a model with 
absorbing and refiecting strips in DW96). 

The dynamics of an irregular grain subjected to 
isotropic radiation is very similar to a grain subjected to 
the Purcell's torques arising, for instance, from H2 forma- 
tion. For instance, one would expect to have thermal trap- 
ping of sufficiently small grains due to thermal fiuctuations 
as described in LD99a. Therefore RATs induced by isotropic 
radiation (henceforth "isotropic RATs") only marginally al- 
ter the problems that the paramagnetic alignment mech- 
anism faces in explaining observational data. In addition, 
as we mentioned above, the "isotropic RATs" are usually 
weaker than those that arise when a grain is subjected to 
anisotropic radiation (see DW97). 

Due to the situation described above, for the rest of the 
paper we shall associate RATs only with the part arising 
from anisotropic radiation, as, for instance, was done in Cho 
& Lazarian (2005). In other words, we treat the torques 
arising from isotropic radiation as a particular realization of 
the Purcell's torques. 

Let us now introduce briefiy some basic definitions. 
RAT Trad is defined by 



where Qr is the RAT efficiency, 7 is the anisotropy degree, 
and u\ is the energy density of radiation field of the wave- 
length A. Here Oe// is the effective size of the grain which 
is defined as the radius of a sphere of the same volume with 
the irregular grain (similar to DW97). In general, Qr is a 
function of angles O, /3, $ in which Q is the angle between 




Scattering coordinate systsm 



Figure 1. The orientation of a grain, described by three principal 
axes 0,1,0,2,0,3, in the laboratory coordinate system (scattering 
reference system) 61,62, 63 is defined by three angles ©, /3, <I>. The 
direction of incident photon beam k is along ei . 

the axis ai corresponding to the maximal moment of in- 
ertia (henceforth maximal inertia axis) with respect to the 
radiation direction k, /3 is the rotation angle of the grain 
around ai, and "1> is the precession angle of ai about k (see 
Fig. [1} . To help the reader familiar with the earlier works 
on RATs, wherever possible, we use the same the notations 
as in DW96 and DW97. 

The RAT efficiency can be decomposed into compo- 
nents in the scattering system via 

Qr(e,/3,<E>) =Qei(e,/3,0)ei 

+ Qe2(e,/3,0)(e2 cos <E>-f 63 siu $) 

+ Q,3(e,/3,0)(e3 cos $-62 sin $), (2) 

where ei, 62, 63 are shown in Fig. (2] In addition, for the sake 
of simplicity, we have denoted Qei(0, P, 0) = Qr(6, P, 0).ei, 
Qe2(e,/3,0) = Qr(e,A0).e2, Qe3(e,/3,0) = 
Qr(6, /3, 0).e3. In what foUowings, we use Qei,Qe2,Qe3 
for the RAT components and keep in mind that they are 
functions of at $ — 0. However, in some particular 
cases, these angles will be explicitly written. 



3 INTRODUCING ANALYTICAL MODEL 
(AMO) OF A HELICAL GRAIN 

Let us consider an asymmetric grain shape consisting of a 
refiecting spheroid and a square mirror with the side h at- 
tached on a pole of the length h ■ For the sake of simplicity, 
we assume that the mirror and the pole are weightless. Also, 
both the mirror and the spheroid are assumed to be perfectly 
refiecting. Moreover, we neglect the shadowing of the mirror 
by the grain by assuming that h ^ I2 (see Fig. 

3.1 RATs from a reflecting spheroidal body 

Consider first RATs acting on an oblate spheroidal body. As 
a consequence of its symmetries, such a grain is not expected 
to exhibit any spin-up arising from RATs, provided that the 
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View fi'om tlie right 





Light beam ^Bf] 


















H k 



= 



Figure 2. Our grain model consists of a mirror connected to 
an oblate spheroid by a weightless rod. The distance between 
the mirror and the spheroid is assumed to be much larger than 
the mirror's size. Both the mirror and the spheroid are perfectly 
reflecting. 



incident radiation is not circularly polarized. The latter will 
be our assumption for the rest of the paper. 

Consider a photon beam of wavelength A, propagating 
in the k-direction parallel to the axis ei of the lab coordinate 
system (see Figs [1] and [2} . The momentum of a photon is 
deposited to the grain as it reflects from the grain surface. 
As a result, the grain experiences a net torque. 

For simplicity, we assume that the grain rotates fast 
around the maximal inertia axis, so averaging over such a 
rotation is suitable. Also, from equation ([2]) it follows that we 
only need to find RATs for $ = 0. Therefore, the resulting 
RAT is 

J^rad=-^—^ {QelGl + Qe2e2 + QeSBs), (3) 

where b is the length of the major axis of the spheroid, RAT 
components are given by (see Appendix A for their deriva- 
tion) 



Qe3 

Qel = 0, 

Qe2 = 0, 



2ea 2 



l)K{e,e) sin 29, 



(4) 

(5) 
(6) 



where s — a/b < 1 is the ratio of the minor to the major 
axes, K[Q,e) is a fitting function depending on O and the 
eccentricity of the oblate spheroid (see the lower panel in 

Fig. Eg. 

Following equations (j4|-([6l we see that a reflecting 
spheroidal grain does not produce any Qei,Qe2, but 
only. Therefore, it is easy to see that the only effect of RATs 
on the spheroidal grain is to cause the precession in the plane 
perpendicular to the radiation direction k||ei. 



3.2 Torques from a reflecting mirror 

Consider now torques that act upon the perfectly reflect- 
ing mirror attached at an angle to the oblate spheroid (see 
Fig. ^ . The pole is considered too thin to interact with the 



radiation0 The normal unit vector N which determines the 
orientation of the mirror in the grain coordinate system is 
given by 



N = niai -I- n2SL2, 



(7) 



where ai, a2, a.-s are the principal axes of the grain (i.e., the 
principal axes of the spheroid because the mirror and the 
pole are weightless). Here ni = sin a, n2 = cos a with a is 
the angle between N and a.2 (see Fig. [SJ. 

Due to the rotation, the cross section of the mirror with 
the surface area A, varies as (see Appendix B) 

A± = yllei.N| = A\ni cos 6-112 sin 6 cos /3|, (8) 

Following the same above procedure (see Appendix B 
for detail), we get RAT 



(Qeiei 4- Qe2e2 + Qeses) 



(9) 



where I2 is the size of the square mirror, h is the length of 
the pole, and RAT components are given by 

Qel = -T-l"-i cos B — 712 sin B cos /3j[nin2 cos^B 
A 

+ — cos (3 sin 2B — cos f3 sin 2B 

— nin2 sin^B cos^/3], (10) 

4^1 2 2 

Qe2 ~ — T-|"-i cos B — 712 siu B COS /9| [n^ cos /3 cos B 
A 

cos p sm 2B — sm 2B + 712 cos p sm BJ , 



2 

ill, 



(11) 



Qe3 = — -t-Itii cos B — 712 siu B COS /3|ni sin f5[ni cos B 
A 

— 712 cos /3 sin B]. (12) 
For B = 0, TT, the RAT components for the mirror are 
4/inf7i2 



A 

4/in? 

4/in? 
A~ 



cos /3, 
sin /3, 



(13) 
(14) 
(15) 



Equations (|13p . (|14|l and (|15[) reveal that Qei does not de- 
pend on l3, but Qe2, Qei are periodic functions of /3. As a 
result, when averaging over the rotation angle /3 from to 
27r, Qe2, Qe3 vauish for B = 0, tt (see a proof based on more 
general symmetry considerations in ^4]2}. 

For arbitrary B, assuming that the rotation of the grain 
around the shortest axis is very fast, we can average RATs 
over P from to 27r. The resulting RAT components are 
given by 

47rZi7iin2 



Qe 



A 



-(Scos^B- l)/(B,a), 



47rZi7iin2 . „^ X 

Qe2 = r sm 2Qg(Q,a) 

A 

Qe3 = 0, 



(16) 

(17) 
(18) 



^ In our model, the only purpose of the existence of the pole is 
to minimize the effects of shadowing of the oblate grain core by 
the mirror. 
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where f{0, a), g{Q,a) are fitting functions depending on 
0,Q wliicli cfiaracterize tfie influence of variation of cross 
section on RATs. Tlie analytical approximations for them 
are given in Appendix B. Note, that the dependence on 
li arises in equations (|16p and l|17|) due to the assumption 
A <C ^1. In the opposite limit we expect A to act as an arm 
for torques and therefore no dependences of the RAT effi- 
ciencies Qe\,Qe2 on h/X to cxist. 



3.3 AMO: RATs' Properties 

For the sake of simplicity, for the rest of the paper, apart 
from the Appendices B2 and B3, we consider AMO for a 
single value of angle a — 7r/4. Combining RATs produced 
by the reflecting oblate spheroid (see equations [JJ |S} , and 
RATs induced by the reflecting mirror (see equations 1161 
lisp , for a = 7r/4 AMO has the following components 



]ei = (3cos e - l)/(e,7r/4), 

47r/inin2 . , . . 

?e2 = ^ sm 265(6,71/4), 



A 

2ea{s^ - 1) 
A 



K{e,e) sin 26, 



(19) 
(20) 
(21) 



where the analytical and numerical fltting functions 
/(6,7r/4) and g(6,7r/4) are shown in Figs Bl and B2, re- 
spectively. 

However, as we see further, Qes does not affect the 
alignment apart from inducing the precession. To roughly 
estimate the latter, one does not need to know the exact 
form of A'(e,e) (see p. 

Equation (|19|l reveals clearly that Qei is symmetric, 
while equation (|20p shows that Qe2 is asymmetric. In ad- 
dition, Qe2 is zero for 6 = 0, vr and it/2, i.e., when the 
maximal inertia axis ai is parallel or perpendicular to the 
light direction (see also Fig. [Sjl . 

In particular, two first components, Qei , Qe2 depend on 
the product of projections of the normal vector, i.e., 711712, 
that suggests us to define the helicity of a grain. A grain 
of right helicity is defined so that, when the maximal iner- 
tia axis ai is parallel to the radiation beam, it can rotate 
in clockwise sense around the radiation beam; a grain has 
left helicity if it rotates anti-clockwise around the radiation 
beam. With these definitions, we see from equations (|19p 
and (|20p that the grain with right helicity corresponds to 
the mirror being oriented so that nin2 > 0, and nin2 < 
for the grain with left helicity. It is straightforward to change 
the grain helicity from right to left by a rotation of the mir- 
ror over an angle of n/2. Thus, grains with left and right 
helicity are mirror symmetric, as expected. Naturally, the 
helicity of a spheroid is zero. 

Fig. O shows components of RAT normalized over the 
maximum of |Qei|, so that they exhibit the functional de- 
pendence of RATs on 6. Since Qe3 arises from the reflecting 
spheroid, not subject to the mirror size, in order for AMO to 
be self-consistent, we normalize Qes so that Qea""" = QTi^"^ ■ 
It is evident from examining the upper and lower panels of 
Fig. 13 that the transition from left-handed to right-handed 
grains induces a simultaneous change of the both Qei and 
Qs2 components. This can be easily understood based on 
the fact that both Q^i and Qe2 depend on the product ni7i2, 
which defines the helicity of grain. Therefore, as the helicity 



AMC, a = 45 , Right helicity : 
Q, 




AMO, a = 45 , Left helicity 

Qe, 

Qe2 

Qe3. e = 0.5 




Figure 3. The components of RAT normalized over the maxi- 
mum of I Qei I as a function of for right handed (Upper Panel) 
and left handed [Lower Panel) models of AMO. The helicity 
change is obtained by twisting the mirror by 7r/2. The maximum 
of Qe3 is also normalized to be equal Q^"^. Figs exhibit zeros of 
Qe2, and symmetries of RATs that we study for irregular grains 
using DDSCAT in g] 



changes, i.e., nin2 reverses the sign, Qei,Qe2 change syn- 
chronically. We shall see that this property is also present for 
arbitrary-chosen irregular grains that we study numerically 
in ^ Fig.[3]also shows clearly that the third component Qes 
that arises from the spheroidal body of the grain does not 
change as the grain gets the opposite helicity. This exactly 
what we expect from a spheroid. 



4 RATS: AMO VERSUS IRREGULAR GRAINS 

We justify AMO's utility finding the correspondence of the 
functional form obtained for the torques that our toy model 
experience with those exerted on actual irregular grains. We 
start with finding generic properties of RATs for irregular 
grains using general symmetry considerations and follow fur- 
ther with numerical calculations. 
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4.1 Symmetry considerations for RATs 

Here we show that some properties of RATs foUow from gen- 
eral considerations based on the analysis of symmetries. For 
instance, we have observed that for AMO Qe2, Qe2 become 
zero after /3-averaging at points = and vr, while Qei 
does not depend on /3. This property is valid for arbitrarily 
shaped grains. Indeed, when Q is either or tt the radia- 
tion direction presents the axis of symmetry. It is obvious, 
therefore that changes of j3 cannot change the RAT com- 
ponent along the radiation direction (i.e. does not change 
Qei), while any perpendicular component of RAT (i.e. both 
Qe2 and Qes), should vanish as the result of /3-averaging 

Further on, we discuss the properties of /3-averaged 
RATs. We can observe, that, similar to the case of AMO, 
the component of RATs Qei for irregular grains is symmetric 
with respect to ^ tt — change (see Fig. [S| . This sym- 
metry is not exact, but it gets better for grains for which 
mutual shadowing of dipoles gets less. The symmetry of Qei 
ensures that the torque along k has the same sign and sim- 
ilar magnitude when the grain flips over. At the same time 
the RAT component Qe2 is anti-symmetric (see Fig. |6} , it 
changes the sign for a transformation tt — 0. This also cor- 
responds to AMO. Similarly to Qe3 the symmetry of Qei is 
only approximate. 



Comparison 





: ' AMO', left helicity 
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Figure 5. Figure shows the comparison of RATs normalized over 
the maximum of |Qei| between AMO (the left handed grain) and 
DDSCAT (for shapes 2 and 4 and monochromatic radiation of 
A = l.2fim). Solid and dashed lines show normalized RATs cor- 
responding to Q^^/Q^"^ = 0.78 and 1 in which the functional 
forms are obtained from the analytical approximation given by 
equations II19I I and I I20II with tabulated functions /, g. Dot and 
dashed-dot lines show normalized RATs for shape 2 and shape 4, 
respectively. 



4.2 Zero points of Qe2 at = 7r/2 

For AMO, according to equation (|21[) . Qe2 is equal to zero 
for — n/2. A similar property also exists for irregular 
grains, as obviously seen in Figs |6] and [T] There it is shown 
that when the maximal inertia axis is perpendicular to the 
radiation direction, the magnitude of Qe2 is very small. This 
can be explained in terms of the interaction of the electric 
dipoles with the electric field vector of radiation as follows: 
the interaction between electric field and electric dipoles in- 
duces their rotation around e2 to emit circularly polarized 
photons. However since ai is perpendicular to k, the elec- 
tric field is only able to induce the rotation of the electric 
dipoles in a plane containing e2. As a result, the torque com- 
ponent vanishes. Quantitatively, according to equation (C5) 
(see Appendix C), it follows that, when — 90", we have 

Qabs.e2 ~ k cos f3Re'^[pj.Einc]e'"^^ . (22) 

Since Xj = rsin asin /3, it follows 

Qabs-e2 ~ ^ cos /3 cos (A sin l3)[pj .Einc], (23) 
j 

where the term [pj.Einc] is a function that is independent 
of /3. It is obvious that Qa6s.e2 is a function of P which 
is zero when averaging is performed for /3 over [0, 27r]. We 
calculated Qaba-e2 for difi'erent A/a and grain shapes, and 
found that Q^2 is indeed close to zero at cos0 = 0, which is 
consistent with our analytical expectation. 



In fact, these considerations prove not only the zero values of 
Qe2 and Qe3 but also their periodicity as a function of f). 



4.3 DDSCAT Calculations 

Fig.|4]presents the test grain shapes that we have calculated 
RATs for using DDSCAT. Parameters for calculations are 
given in Table[T] Shapes 1, 2 and 3 have been used in DW97. 
We added to them shapes 4 and 5. In addition, we created 
a mirror symmetric shape of shape 1, namely, shape 1* and 
provided the DDSCAT calculations for a spheroidal grain 
(see more details in Table [!}. We adopt dielectric functions 
for astronomical silicate in which a feature in the ultraviolet 
is removed (see DW97; Weingartner & Draine 2001; Cho & 
Lazarian 2005). 

We discussed for AMO, that the sign of helicity can 
be changed by taking the mirror image of the grain. We 
performed a similar procedure to the irregular grains and 
obtained results similar to the ones obtained for AMO (see 
Figs. [7] and [8]). 

Note, that we observe that Qei and Qe2 change syn- 
chronously when we calculate torques for a mirror image of 
a grain. We see that the shape 1 has one type of helicity, 
while shapes 4, 5 and mirror symmetric image of shape 1, 
i.e. shape 1*, have another type of helicity. 

Fig. [5] provides a comparison of normalized RATs be- 
tween AMO and DDSCAT calculations performed for two 
irregular grains induced by monochromatic radiation field of 
A = 0.2^m (see more in i]4.4p . It can be seen that they pos- 
sess the same symmetric properties as well as zero points. 
Also, the functional form of normalized Qei and Qe2 cal- 
culated for the irregular grains and AMO are remarkably 
similar, in particular for the Qe2 component. Typically, the 
RAT components for shape 2 are similar to those of AMO 
with Q^-C'' IQ^r = 1 ratio, while RATs of shape 4 are sim- 
ilar to those of AMO with Q™'^/QS"'' = 0.78 . Hence, by 
changing the ratio of amplitudes of the RAT components 
for AMO, we can obtain analytical expressions of RATs for 
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Figure 4. Geometry of grains under study: shape 1, 2, 3 are similar to those of DW97, shape 4, and 5 are created from 15 and 11 cubic 
blocks respectively, and an ellipsoidal shape. 

Table 1. Grain shapes and parameters for calculation of RATs 



Grain shapes Dipole # Size (fim) Wavelength {fim) Helicity 



Shape 1 


832000 


0.05 - 0.2 


ISRF 


right 


Shape 1* 


53248 


0.2 


1.2 


left 


Shape 2 


45056 


0.2 


ISRF 


left 


Shape 3 


102570 


0.2 


ISB.F 


left 


Shape 4 


15000 


0.2 


1.2 


left 


Shape 5 


11000 


0.2 


1.2 


left 


Hollow 1 


832000 


1.0 


0.1 


right 



For all calculations here, we adopt the dielectric function for astronomical silicate. 



a number of irregular grains. To have RATs appropriate to 
irregular grains, it is necessary to use DDSCAT to estimate 
the magnitude of RATs. Combining functional forms from 
AMO and magnitude from DDSCAT, we can obtain an- 
alytical approximate expressions for RATs components of 
irregular grains. Note, that in Figs [3] and (5] we normalized 
RATs over IQ™""! = IQei(e = 0)], that gives rise to Q™i"^ 
remained the same for all realizations of AMO and the ir- 
regular grains. It is easy to see that with this choice, AMO 
reproduces very well Qe2 for irregular grains, but gets Qei, 
which is a bit larger at cos O = (O = 7r/2) than that 



for irregular grains. Potentially, this may mean that more 
appropriate parametrization should include Q™"^, which is 
not defined as Qei at | cos 0| = 1 (O = or tt) as we do 
in this paper, but, for instance, the amplitude value of Qei, 
which is |Qei(e = 0)| + lQei(e = 7r/2)|. We feel, however, 
that the our present parametrization has the advantage of 
simplicity and is sufficiently accurate. 
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Figure 6. Qel and Qe2 for grains of left ( Upper Panel) and right 
{Lower Panel) helicity. The symmetry of Qel with the transfor- 
mation TT — is clearly seen for all grain shapes. Qe2 is antisym- 
metric with respect to the same transformation tt— 0. A similarity 
with the torques produced by AMO is evident (see also Fig|5j- 

4.4 Parameter study 

Above we compared the properties of RATs in AMO with 
those obtained numerically from DDSCAT for a few chosen 
grain shapes and radiation spectra. To see how general our 
numerical results are, we attempt a limited parameter study, 
namely, we study how the properties of RATs vary with the 
spectrum of the incident radiation for different grain shapes. 
One can view the AMO formulae as a physically motivated 
fit to RATs acting on astrophysical grain with Usff/X < 1. 
The parameter study is intended to find out how good is 
this fit. 

Fig. [7] shows Qel and Qe2 for the shape 1 produced by 
different radiation fields. There the upper panel show that 
when monochromatic radiation fields of A/oe// increases, 
the symmetry of Qei and zeros of Qe2 do not change. How- 
ever, their amplitude decreases. In addition, the symmetric 
property of Qei and zeros of Qe2 also remain unchanged 
when being averaged over different radiation spectrum (sec 
the lower panel in Fig. [7]). 

Now let consider RAT properties for different irregular 
grains. The upper panel in Fig. |8] shows RATs for different 



Figure 7. RATs for the shape 1 corresponding to various X/a^ff 
( Upper Panel) and RATs averaged over a range of wavelengths 
{Lower Panel). 



shapes: shape 1* which is a mirror symmetric copy of shape 
1, and shapes 4 and 5 are built from 15 and 11 cubic blocks, 
respectively. We also see clearly that Qei exhibits the sym- 
metry, and Qe2 exhibits the asymmetry that we have already 
seen with AMO and other grain shapes. 

From Figs [7] and [S] it follows that the form of Qei, Qe2 
for shape 1 is mirror-symmetric to the corresponding RAT 
components applied to shape 1*. This mirror symmetry is 
also evident when we compare Qei , Qe2 with those of shapes 
4, 5 (see Fig. |4]). This implies, similar to AMO, irregular 
grains may be of right and left helicities. A comparison be- 
tween Figs |6] [7]iipper, and IHltpper shows that shapes 1*, 2, 3, 
4, 5 are of left helicity, while the shape 1 is of right helicity. 
Note, that Fig. [Tfoiuer clearly shows that the helicity is in- 
dependent of wavelength and is intrinsic attribute of grains 
which is associated to their shape. Especially, we can obtain 
a grain with the opposite helicity by performing a mirror 
symmetric transformation, which is illustrated by AMO in 
Fig. [3] We remind the reader, that the correspondence be- 
tween Qel and Qe2 for AMO and shape 2 and 4 is illustrated 
by Fig. m 

For the third RAT component Qes , it exhibits analogous 
properties with Qe2 obviously seen in the lower panel in Fig. 
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Figure 8. Upper panel: RATs for different grain shapes. Shape 
1* corresponds to mirror-symmetric image of shape 1. Symmetric 
features of Q^i and zeros of Qe2 are clearly found. Lower panel: 
The third component of RATs Qe3 for different shapes are shown 
together with that of a spheroid with e = 0.5 predicted by AMO. 
An analogy exists between the zeros of Qe2 and Qe3- Also, the 
shape of Qe3 is not affected by the change of grain helicity. 



[H Also, for an axisymmetric shape, i.e., spheroid, Qe3 is 
still significant. Note, that its functional form obtained by 
DDSCAT is very consistent with that predicted by AMO in 
see 33] (see dot and dot-dashed line in Fig. [8} . Furthermore, 
we see that it has similar forms for left (shape 2 and 3) and 
right (shape 1) helical grains (see Fig. ^ower). While the 
dependences of Q^i and Qe2 undergo a transformation when 
shape 1 is substituted by the shape 1*, having the opposite 
helicity, a comparison of lower and upper panels of Fig. [S] 
shows that the shape of Qe3 component stays the same. This 
is expected, as Qes does not depend on the helicity of grains. 

Also, the Qe3 component has zeros at = 0, tt, similar 
to Qe2 (see the lower panel in Fig. However, the anti- 
symmetry of Qe3 is less prominent than for Qe2- 

As we will see that the ratio QS'"'/Q3'''^ is an impor- 
tant parameter that determines the existence of high-J at- 
tractor points. For irregular grains, this ratio is a function of 
the ratio of wavelength to grain size, as shown in Fig. [5] for 
three irregular grains. The peak of QS" 



10.0 



E 



1 .0 




0.1 



shape 1 , Ogff=0.2 jJ,rri 
shape 2, o^^<=0.2 /xm 
shape 3, agff=0.2 fj.m 
shape 1, agff=0.05 /im 
shape 1 , ag,f = 0. 1 /j.m 



10 



A/a. 



Figure 9. Ratio of Q™'^,Q™^ as function of wavelength to 
grain size X/a^ff for different grain shapes and sizes. The shaded 
region corresponds to the parameter space in which the high-J 
attractor point is present for the alignment in respect to the beam 
direction (see ^El . 



for different shapes. The form of the curve for sufficiently 
large ratios of X/a^ff can be approximated as 10 ""^^ . This 
dependence can be used to reduce the number of DDSCAT 
calculations necessary for determining the alignment for ar- 
bitrary radiation fields. 

In general, we observe strong similarities between the 
plots of RAT components obtained for very different grains. 
Thus we can expect that the RAT alignment should be sim- 
ilar for such grains. We will discuss the alignment for AMO 
and irregular grains in § [6] and § [T] 



4.5 RATs: Comparison with AMO 

The actual grains are not perfectly reflecting particles and 
the scattering that they induce cannot be described by geo- 
metric optics that we employ for AMO. Therefore the justifi- 
cation of the AMO utility can be obtained via a comparison 
of the functional form of the torque components obtained for 
irregular grains with the AMO predictions. Naturally, one 
should not expect to see the amplitudes of the torques to be 
the same. Therefore the comparison should be done for the 
normalized torque components. However, we preserve the 
ratio of the components. 

Naturally, our sample of RATs acting on grains stud- 
ied with DDSCAT is limited. It includes several grain sizes. 
For instance, for shape 1, we studied for grain sizes of 
0.05,0.08,0.1 and 0.2/im. For other shapes, the size fle// = 
0.2/im is studied, except the hollow shape 1 with fle// = 
l^m. We calculated RATs for the entire spectrum of ISRF 
corresponding to 21 wavelengths in the range A = O.l/im to 
lOO^m, for shape 1, 2 and 3, and the monochromatic radia- 
tion with A = 1.2^m for other shapes. This provides us with 
RATs calculated for 130 realizations of grains and radiation 
fields Q. This makes it the most extensive sample of RATs 



* Each realization corresponds to a given grain size and a given 
wavelength 
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studied numerically. It is obvious, that in our paper we can- 
not present plots of the RATs for all the realizations that 
we calculated (see e.g. Fig. [5)|. A quantitative comparison 
based on the deviation testing for normalized RATs of all 
realizations and AMO will be presented below. 

We show in Fig.[9]that, for irregular grains, the relative 
amplitude of Qei versus Qe2 changes both with the grain 
shape and wavelengths. However, our studies in this paper 
shows that the functional form of the RAT components for 
all the cases we studied is still well represented by AMO 
(with different ratio of QTi""" / , e.g.. Fig. [g. In other 
words, while DDSCAT studies of alignment for grains of a 
few chosen shapes cannot reveal the generic properties of the 
RAT alignment, revealing the correspondence of the func- 
tional dependences of the torques between irregular grains 
and AMO provides a deep insight into the alignment. 

Since we are only interested in the functional forms of 
RATs, let us introduce the mean deviation over O for the 
components Qei and Qe2 as foUowings 



{A')(Oe 



(Q 



DDSCAT 



(6) 



)^."°(e))^de, 

(24) 



where 



(9) denote Qei(O), (3e2(0) for irregular 
grains, Q^i^'^iQ) denotes the torque components for AMO 
in which the relative magnitude are rescaled to have the 
same ratio Q^-^"^ IQ^^^ with each realization of irregular 
grainijfl- In equation (|24p . Q^i"'^ is the maximum of Qei and 
(3e2, that is chosen the same for both AMO and irregular 
grains. 

We perform (A^) testing for our sample consisting of 
130 realizations of irregular grain shape, size and wave- 
length. To see the correspondence of AMO with different 
grain shape, size and wavelength, in Fig. [TU]we show (A^) 
as a function of A/fle//- 

Fig. llOb pper shows a good correspondence for the com- 
ponent Qei between irregular grains and AMO. The value of 
(A^) ranges from as small as 10"'' to 2 x 10~^. In addition, 
Fig. llOt oioer shows an extremely good correspondence for 
the component Qe2 between all cases of irregular grains and 
AMO. The value of for different sizes of shape 1 are nearly 
the same, but it changes with grain shapes. For instance, 
shape 2 and 3 have the better correspondence with AMO 
than shape 1. The corresponding value of (A^) is about 
2 X 10^^ and 3 x 10~^ for shape 2, 3 and 1, respectively 
(see Fig llOH ower'). For both components and all grain shapes 
and sizes, (A^) do not change so much with respect to wave- 
length of radiation field for A/ae// > 3. However, it increases 
when A/fle// decreases. It is clear that the smaller (A^) is, 
the better correspondence is. Therefore, the observed worse 
correspondence for the range A/fle// < 3 can been easily ex- 
plained in terms of the RAT properties for irregular grains. 
In fact, following Fig. [7]upper, it can be seen that the com- 
ponents Qei, Qe2 for A/oe// = 2.8 are indeed less symmetric 
than for A/fle// = 8 and 15. In other words, the symmetry 
of RATs decreases toward small A/ae//. 



^ We note again that, throughout this paper, apart from Appen- 
dices B2 and B3, the functional forms of RAT components for 



AMO corresponds AMO with a 
maximum is adjustable. 



45*^, and the ratio of their 
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Figure 10. Figs show (A^) testing as function of X/a^ff for 
normalized Qei (upper panel) and Qe2 (lower panel) between 
irregular grain shapes and sizes with AMO. 



Note, that lower values of (A^) can be achieved for small 
A/se// ratios, if Qei is considered separately for O within 
ranges [7r,7r/2] and [7r/2, 0]. This procedure can compensate 
for the differences in the surface characteristics on the side 
of the grain towards the light direction and opposite to it. As 
we discuss in ^such differences are natural for less idealized 
models of helical grains. We do not persue this approach 
within this paper, prefering a simple model with a reasonable 
fit to a more complex model that can provide a better fit. 
Nevertheless, in future, dealing with particular cases, e.g. 
graphite grains, the approach of considering different ranges 
of Q may prove advantageous. 

Therefore, our quantitative comparison of the RATs 
based on (A^) for irregular grains with AMO implies that 
the functional form of RATs obtained in the limit A ^ ae/ / 
(AMO) is also valid for RATs in the limit A ~ Ce// and 
A > ae// (for irregular grains). This allows us to use AMO 
as a representative model for describing RATs of realistic as- 
trophysical grains, e.g. study grain alignment using AMO. 
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5 AMO: DEFINITION AND 
GENERALIZATIONS 

Our studies above show that the model of AMO correspond- 
ing to a = 7r/4 corresponds well to the results of numerical 
calculations of RATs, provided that we treat Q^^^'/QS"^ 
as an adjustable parameter. Indeed, while we show in Ap- 
pendix B that this ratio changes with a, the range of the 
ratio variations is lower than for irregular grains (compare 
Fig.|B6]and FiglJ. 

This adjustment is natural, as we do not really expect 
the scattering by irregular dielectric grains to be entirely 
equivalent to scattering by mirrors of our toy model of a 
grain. 

Our analytical studies are based on AMO, which is pro- 
vides a simple model of RATs. If necessary, AMO can be 
trivially generalized by adding additional mirrors to grain 
surface and accounting for a partial shadowing of the mirrors 
by the spheroidal grain body. In addition, one can consider 
not perfectly reflective models, but models with refractive 
grain body and also refractive plate instead of a mirror. This 
allows us to vary AMO's properties. 

An arbitrary attachment of mirrors may make Qei less 
symmetric. For instance, a grain shown in the upper panel 
of Fig. [11] has a mirror attached to one of its surfaces. Nat- 
urally, turning this surface towards the beam of radiation 
produces torques that are different from the case when the 
mirror is hidden by the elliptical body of the grain. 

In addition, the adopted AMO does not demonstrate ac- 
count for torques arising from absorption. If, however, one 
adds to the spheroidal grain a one-dimensional damped os- 
cillator with the axis at an angle to the spheroid symmetry 
axis (see the lower panel in Fig. Ill[) . such a grain would 
have a non-zero torque component Qe2 arising from light 
absorption. Note, that just adding absorbing component to 
one of the surfaces of the reflecting mirror would distort the 
symmetry of torques Qei, which would correspond to RATs 
calculated for irregular grains for small X/a^ff ratio. 

While adding such effects could provide us with better 
fits of RATs obtained numerically for particular shapes, this 
would increase the complexity of our model. We do not pur- 
sue this path here, therefore. However, we can see that in 
some cases, e.g. dealing with highly absorbing grains, e.g. 
graphite grains, one may have to make the corresponding 
analytical model more sophisticated. 

Note, that some modifications do not require making 
the model more complex. Historically, spheroidal grains were 
used to demonstrate alignment. However, for demonstrat- 
ing the effects of H2 torques Purcell (1979) considered a 
"brick-shaped" grain. If we want to study the effects of H2 
torques for AMO, it is natural to consider not a "spheroid 
with a wing", as we have done in this paper, but a "brick 
with a wing". In terms of AMO this would not change the 
properties, apart from the value of Qes torque, which is not 
important for the alignment, anyhow. 

Intentionally, our AMO is based on an intuitive macro- 
scopic model. Microscopic models based on the analytical 
treatment of RATs are possible also. For instance, light 
scattering by optically active (chiral) sphere was analyti- 
cally studied by Bohren (1974). The author derived ana- 
lytical solutions for scattering matrix. Due to the optically 
activeness, such spherical grain can produce RATs. This 



Spheroid + a shadowed mirror 




A damped oscillator and a spheroid 




Figure 11. Models of a grain shape consisting of a mirror shad- 
owed partially by the oblate spheroid (upper panel) and a par- 
tially damped oscillator attached to the oblate spheroid (lower 
panel). 



model would correspond to the Dolginov (1972) suggestion 
of quartz grains being spin up and aligned by RATs. 

An important generalization of AMO is to consider that 
its inertial properties are given not by a spheroid but by a 
triaxial ellipsoid. We expect changes to result in difference 
in dynamics during periods of crossovers. The corresponding 
effects together with the modification of inertial properties 
of AMO are considered in Hoang & Lazarian (2007). 

All in all, while it is rather easy to make AMO more 
sophisticated and provide a better fit for RATs within dif- 
ferent ranges of A and a^ff, for the rest of the paper, we 
adopt a very simple model of AMO. This model provides a 
reasonable flt to generic properties of RATs acting on actual 
irregular grains. Thus studying grain alignment with AMO 
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should provide insight into the alignment processes acting 
upon astrophysical grains. 



6 ALIGNMENT WITH RESPECT TO LIGHT 
DIRECTION 

While the earlier studies dealt with the alignment in respect 
to magnetic field, in this section we show that RATs can 
align grains on their own, without any influence of magnetic 
fields. In this case, the direction of radiation k is the axis 
of alignment. As we further discuss in ^9.11 such alignment 
happens in the presence of magnetic field when the rate of 
precession induced by radiation is faster than the Larmor 
precession rate. Dealing with this simple case also prepares 
us for dealing with a more complicated cases of alignment 
in respect to magnetic field in ^ 

6.1 RATs: spin-up, alignment and precession 

To understand the role of the RAT components, we calculate 
torques that spin up, align, and induce grain precession. A 
RAT component that acts to spin up grains, H, is directed 
along ai, the component that aligns grains, F, is directed 
along Q, and RAT that causes precession, G, is along <t. 
They are respectively given by (see DW97) 

H{e, $) = Qp.ei(e, $)cose + Qr.e2(e, $)sinecos<l> 

+ Qr.e3(e, $)sinesin$, (25) 

F{e,^) = -Qp.ei(e,$)sine + Qp.e2(e,<l')cosecoscE> 

+ Qr.e3(e,<l>)cosesin$, (26) 

G(e, <!>) = -Qr.e2(e, $)sin$ + Qr.e3(e, $)cos$. (27) 

On the other hand, following equation ([2]), RATs at a 
precession angle "1> are related to which at $ = via 

Qr.ei(e,$) = Qei(e,0), (28) 
Qr.e2(e,$) = Qe2(e,0)cos$-(5e3(e,0)sin$, (29) 
Qr.e3(e, $) = Qe2(e, 0)sin$ + Qe3(e, O)cos$. (30) 

Plugging the above equations into l|25[l . I|26p and (|27p 

we get 

G(e,$) = Q.3(e,o), (31) 

H{Q,$) = Qei(e,o) cose + Qe2(e,o) sine, (32) 

F(e,-i') = -Oei(e,o) sine + 0e2(e,o) cose. (33) 

If for G(e, $) we are mostly interested in its amplitude, 
the functional form of F{0, $) and H{0, <E>) is essential for 
grain alignment. The problem is that F{&,$) and H{Q,$) 
as well as their counterparts obtained in the presence of mag- 
netic field (see equations [68] and [69]), vary substantially from 
one grain to another. For AMO different grains correspond 
to different ratio As we mentioned in ^4.31 and 

MA\ irregular grains are different in terms of RATs for the 
radiation of different wavelengths and different grain sizes. 
However, for both AMO and irregular grains, the generic 
properties of the RAT components (i.e., symmetry of Qei, 
as well as the asymmetry and zeroes of Qe2,Qe3) is always 
unchanged. Therefore, unlike Qei and Qe2, the components 
F(e,$) and H{0,$) do not demonstrate a universal be- 
havior and play an auxiliary role in our study. 

We see that the precessing torque depends only on the 



third component Qs3{O,0), while the aligning and spin- 
ning torques are related to two first components, namely 
(5ei(e, 0), (5e2(e, O). We note that the functions F, G, H are 
the functions of only variable O, in this case. 

For AMO with a = 45'^, substituting analytical expres- 
sions / = /^/2 and g given by equations HB28|I and HB31|I 
into equations ([T6|l and (flTj) . we get 



ir^ n\ 16;inin2|n2|, r,N 
l(Q.O) = ^ — ^(5 cos e-2), 



3A 



(34) 



3A 



Qe2(e,0) = 



Therefore, equations (I3HI - (|33| I becomes 

G(e) = -Hff:^!— ilA-(e,e) sin 26, 
A 

, , 8/inin2jn2i , , c m -2^ 

//(O) = — - cos B[l -f 6.91 sm O 

oA 

+ cos^e(5-|- 1.382 sin^e)], 

, , 8/inin2|n2i . , , ^ ni 2^ 

t (B) = — — ^ — - sm e[— 1 + 6.91 cos 6 

oA 



sin 26(1.191-^0.1382 cos^6). 

(35) 



(36) 



+ 1.382cos''e + 5 sin^6]. 



6.2 Simplified Treatment of Crossovers 



(37) 
(38) 



In general, the maximal inertia axis ai of our model grain 
can precess about the vector of the angular momentum J. 
In the present paper, however, for the sake of simplicity, we 
assume a perfect internal alignment, i.e. J||ai. This assump- 
tion coincides with that in DW97 and can be justified by 
the high efficiency of the internal relaxation within a wob- 
bling grain. This relaxation stems from the Barnett relax- 
ation discovered by Purcell (1979) and/or nuclear relaxation 
introduced in LD99b. However, these relaxation processes 
provide a good coupling only when J ^ Jth ~ (kTdli)^^^ 
where Td is the dust temperature, Ji is the maximal moment 
of inertia of grain (Lazarian 1994), i.e. when a grain rotates 
with suprathermal velocities. This condition is not satisfied 
as a grain approaches crossover points, i.e. as J — > Jth- 

We adopt below a simplified treatment of crossovers, 
which is different, however, from the treatment of crossovers 
in DW97. There it was assumed that J||ai up to the mo- 
ment of the angular velocity getting zero. After that DW97 
assumed that J changes its direction to the opposite, while 
the direction of ai is preserved. Such a model of crossovers 
differs from the earlier work on crossovers, which suggests 
that, as the angular velocity goes to zero, the grain un- 
dergoes a flip, i.e., J preserves its direction, while the di- 
rection of ai changes to the opposite (Spitzer & McGlynn 
1979; Lazarian & Draine 1997). Therefore, in what follows, 
we adopt a model in which J||ai up to a crossover; at the 
crossover ai flips, while the direction of J is preserved. This 
makes the dynamics of grains very different from that in 
DW9fl Ahhough we accept that our model is not precise 



^ We may mention parenthetically, that the dynamics that we 
get has similarities to the dynamics of grains with thermal fluc- 
tuations accounted for, as in WD03. There, however, it was in- 
terpreted as a new eff'ect related to thermal fluctuations. 
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Table 2. Physical parameters for diffuse ISM 



Appendix C) 



Definitions 



Typical Values 



Gas density 
Gas temperature 
Gas damping time 
Dust temperature 
Anisotropy degree 
Grain size 
Mean wavelength 
Mean density of ISRF 



n = 30 cm ^ 



: 4.6 X 10l2( 



:J72")^-5 s 

^ gas 

Td = 20 K 
7 = 0.1 
S'cff = 0.2/xm 
A = 1.2 /im 



The normalized values that wc also use are Tgas = 
Tgas/WO K, n = n/30 g cm~'^, and a_5 = a^f f /1Q~^ cm, and 
normalized grain density p = p/3 g cm~^. 

when J ~ Jth , we claim the our simplified model is "roughly 
true" . The latter point is justified for regular crossovers (i.e., 
crossovers without thermal fluctuations) in 381 



6.3 Stationary and singular points 

6. 3. 1 Equations of motion 

To find out whether grains can be aligned with respect to k, 
we need to find stationary states for grain motion exerted 
by RATs. We start with equations of motion 



dJ_ 

It 

— = M ^ 
dt J 



MH((d) - J, 
F{Q) 



(39) 
(40) 



where time t and angular momentum J are scaled in the 
units of the gas damping time tgas (see Table [2]) and the 
thermal angular momentum I\uit- Here I\ is the maximal 
moment of inertia and lot ~ 

[hlp^f/^ is the thermal angu- 
lar velocity with Tgaa is the gas temperature. In equations 
(I39l) and (HOl) 



(41) 



M = 



2I1LUT 



contains important parameters of the physical problem. 
There Urad is the energy density of radiation field, A is the 
wavelength. In the present paper, for our estimates, we use 
the values provided in Table[2]for interstellar medium. It will 
be explicitly stated in case other values of these parameters 
are used. 



dF{e) 



H(Q)dQ 

dF{e) 



< 1, 



< Fie) 



dH{e) 



dO 



(42) 



(43) 



At the stationary points, F{Os 
and (1431) reduce to 



< 0, for HiQs) > 0, 



dF{e) 



dO 

and 
dF{e) 



dO 



> 0, for HiQs) < 0. 



8=03 

0, so equations (|42|l 
(44) 



(45) 



Consider that initially a grain has the maximal inertia 
axis ai parallel to k, i.e., Qo = 0, then F{Oo) = 0, H{0o) = 
Qei(9o)- These assumptions are usually used to estimate the 
grain angular velocity induced by RATs (see DW97; Cho & 
Lazarian 2005). Therefore, the equations of motion become 



at 



J. 



Solutions for the above equations are easily found, 
6 = 0, 

J = MQeiie) + (Jo - A/Qei(eo))e"', 



(46) 
(47) 

(48) 
(49) 

where Jo is the initial angular momentum of grains. Equa- 
tion ((49| shows that as t > 1, J = MQci{0 = 0). In other 
words, grains initially parallel to k are spun up by RATs and 
aligned with k at angular momentum Js = MQei(6s = 0), 
regardless of its initial angular momentum Jo. However, the 
particular point of phase trajectory map may not provide 
stable orientation of the grain, i.e., may not be attractor 
points. [3 

For sin Oo 7^ , substituting F{0),H{0) from equa- 
tions (|37p and (|38|) into equations p9[) and ((40]), we obtain 



(50) 

(51) 

(52) 
(53) 

(54) 



dJ 
'dt 



8Zinin2|n2j 
3A 



M cos e[l + 6.91 sin e 



+ cos^e(5 + 1.382 sin^e) 
dQ _ 8linin2\n2\M 
'dt ~ 



J, 



sin e[-l + 6.91 cos^e 



3AJ 

+ 1.382cos''e + 5 sin^e]. 
Hence, we can easily find a stationary poinl|!| 
=0, 

_ 48Zinin2|n2|M 
•^'^ 3A 
From equations (|37|) and p8|) we get 



dF{e) 



dO 



90.232/1 rtin2|n2 1 
3A 



6.3.2 Stationary points 

The stationary point is determined by setting the equations 
of motion (|39[) and (|40p equal to zero. This point may be an 
attractor or a repellor point. We remind the reader that for 
a stationary point to be an attractor point, it requires (see 



^ The shaded area with diagonal lines in Fig 1241 corresponds to 
the existence of such points. We see that many wavelengths the 
condition is not satisfied. 

* Another stationary point © = 7r, J < is forbidden since J can 
not be negative. Note, that the analysis in this section is applied 
for the right helical grain, i.e., rain2 > 0. 
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and H{Q3 = 0) > 0. As a result, the stationary point Os = 
is a repellor point. 

In a general case, F{0), H{0) are given by equations 
l|32p and (|33p . Thus, substituting them into equations (|44p 
and (|45[) . we get a criteria for a stationary point Os being 
an attractor point, namely, 

-Q'el sin Qs + Q'e2 COS Qs ^ ^ .gg. 
Qel COS 6s + Qe2 Siu 6s 

where Q'^i = ^^^,Qe2 = However, our approximate 

formulae above catches the phenomenon correctly and we 
still get repellor points for 6s = with the exact expressions 
for AMO. 

As the range for ratios of Qei /Qe2 for an arbitrarily cho- 
sen irregular grain may be different from AMO, for some ir- 
regular grains we still may have attractor points for 6s ~ 0, 
as it is shown in Fig [24] for shape 1 subjected to interstel- 
lar radiation field. There the relevant data corresponds to 
tp = 0, as in this situation the direction of magnetic field 
and the light coincide and thus our considerations about 
the alignment in respect to k are applicable. 

6.3.3 Singular points: crossovers 

When J = 0, we have a singular point (see equation I40p . 
In terms of our simplified equations it corresponds to a 
crossover. As we discussed earlier, a crossover is a period 
in which a grain spins down to the point that the compo- 
nent of angular momentum parallel to ai gets zero. In the 
presence of strong internal relaxation that tends to align J 
and ai, this means that J should get small (Spitzer & McG- 
lynn 1979) . Our equations above are derived in the assump- 
tion of J||ai and therefore can not treat grain crossovers 
(cf. §6 [5]). However, they can still trace the grain dynamics 
as the grain phase trajectory approaches the crossover and 
J — *■ Jth- Assuming that initial angular momentum of the 
grain Jo ^ Jth, we disregard the difference between J — Jth 
and J — 0. Thus, from equation (|40p it follows that, in or- 
der to have a physical crossover, grains must have J = 0, 
and the aligning torque F{0) must be zero. The latter con- 
dition is naturally satisfied since we found that F{0) — 
at stationary points. The former one is satisfied as RATs 
act to decelerate the grain rotation. Indeed, equation (|38p 
shows that F{Q) > 0, i.e., it acts to increase 6, for every 
angles 6 < tt. Yet, equation (|37p shows that H{Q) < 
for cos 6 < and H{Q) > for cos 6 > 0. Therefore, 
if initially the angular momentum of grains makes an an- 
gle cos 6 < 0, then their angular momentum is decreased 
due to H{0) < 0; so grains approach directly to the state 
of J = (see the upper panel in Fig. [13] for right helical 
grain). On the other hand, for grains which initially have 
cos 6 > 0, the aligning torque acts to increase 6, while the 
spinning torque increases their angular momentum. Eventu- 
ally, grains attain the angle cos 6 < for which the spin- 
ning torque changes the sign, so they are decelerated to the 
state of J = (see Fig. [13] for phase maps of left and right 
helical grains). 

One can observe that when grains get to the crossover 
with very low J, their maximal inertia axis ai flips with re- 
spect to J to enter the opposite flipping state, i.e., grains 
flip from the upper to lower panel in Fig. \13h mer. for in- 
stance. Right after that, grains flip back to the initial state 



(upper panel). This back and forth flipping process takes 
place frequently. As a result, the crossover point, in our ap- 
proach, can be treated as the attractor point at zero angular 
momentum, hereafter, called low-J attractor poinfl to be 
distinguished from high-J attractor points. 



6.4 RATs: alignment by one component 

The grain alignment for the case at hand is uniquely related 
to two components Qei and Qe2. To understand which com- 
ponent causes the alignment, let us study the role of Qei 
and Qe2 separately. 

If Qe2 = 0, then %^(6s = 0, tt) = 0, equation ([55]) for 
attractor points becomes 

gei(6s = 0)>0, (56) 
Qei(6s = tt) < 0. (57) 

According to AMO, Qei (6s = 0,7r) = 
i6iinin2|n2| j-g^Pg2Q^ - 2) = 1 > 0. As a resuh, we 
expect that the stationary point 6s = is an attractor 
point, i.e., grains are perfectly aligned with respect to k. 
While, the stationary point 6s = tt is a repellor point. 

For the case Qei = 0, besides stationary points 6s ~ 
0, TT as discussed above, there is another stationary point 
corresponding to 6s = 7r/2 (see equationl33p. For 6s = 0, tt, 
equation (|55p reduces to 

dQe 



d6 

dQe2 

d6 



(6s = 0) < 0, 
(6s = tt) > 0. 



(58) 
(59) 



For AMO with a = 45°, thus, 
40hn,n2\n2\ ^-^ 26(1.191 + 0.1382 cos26), 
dQ^IP,_ _ n ^^ - i06.2iinin2|n2| 26s > 0. It mcaus 



Qe2 = 

we have 

<iQc2 tip, _ n •7^^ — io6-2'iiiT 

~de~v^» — U,7rj — 3^ 
that the condition (|58p is not satisfied, but equation (|59p 
is fulfilled. In other words, the stationary point 6s = is 
a repellor points, while the stationary point 6s = tt is an 
attractor point. 

For a particular stationary point 6s = 7r/2, equation 
(I58II in this case becomes 



= < 1. 



(60) 



Q'e2 COS 6s 

Qe2 sin 6s 

Therefore, the stationary point 6s ~ vr/2 is, indeed, an 
attractor point. Note, equation H60p which indicates that 
6s = vr/2 is the attractor point is valid for an arbitrary 
form of Qe2. These results show that the component Qe2 
acts to align grains in the direction perpendicular to the 
radiation beam for arbitrary grain shapes. 



6.5 Phase trajectories 

To show our predictions from the above analysis, we con- 
struct phase trajectory maps in which each RAT component 
acts separately. The difference of our phase maps from those 
in DW97 is that here we present the alignment in respect 



^ In Hoang & Lazarian (2007), we will show that zero angular 
momentum attractor points become attractor points at thermal 
angular momentum due to thermal wobbling 
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to k rather than to magnetic field. In addition, we treat 
crossovers diflerently from how they are treated in DW97. 
With the exception of Figs l23l and l26l to avoid over-crowding 
of our phase maps we do not draw arrows with tgas intervals. 

Throughout the present paper, the phase trajectory 
map represents the evolution of J and angle between J and k 
or Jiifj Each phase map has upper and lower panels, which 
correspond to J parallel and anti-parallel to ai , respectively. 
Also, a circle denotes an attractor point and a cross denotes 
a repellor point. 

Moreover, the title of the trajectory map presents in- 
formation about the model (AMO) or grain shape, size for 
irregular grains. For trajectory maps labeled with AMO on 
their titles, it is a default that the RAT components are 
taken from the exact calculations of equations (|10)I12|I for 
a = 45° which have Q^i'^'/QS'''' = 1-2. For the maps of 
AMO in which the ratio Q'Ti is explicitly shown, it 

means that the functional forms of RATs are similar for the 
case a = 45°, but the relative amplitudes are rescaled to the 
shown value. 

The map for the case when only Qei is present has 
three attractor points: A, B correspond to Qa = 0, and E 
corresponds to zero of Qei at cosO = —0.6 (see the upper 
panel in Fig. I12|l . Therefore, Qei acts to align grains with 
J parallel to k, i.e, ai parallel or anti parallel to k . In 
addition, there are two repellor points C, D. These features 
show clearly our predictions obtained with AMO. 

When Qei — 0, Fig. [12] (lower panel) shows that the 
phase trajectory map has one attractor point D at Oa = n/2, 
which corresponds to a perpendicular alignment, one low-J 
attractor point C at Os = tt, and two repellor points A, B 
as predicted by our above analysis. So Qe2 can align some 
grains with the maximal inertia axis ai perpendicular to the 
light direction. However, we discussed in §2 that Qe2 is equal 
zero at the attractor point Os ~ 7!"/2, while Qei is different 
to zero at this point; thus the alignment with the long axes 
parallel to the direction of light does not occur^3 for the 
case of alignment with respect to k (see § I7.5|l . 

6.6 RATs: alignment by joint action of torques, 
role of QTi''yQT2'"" ratio 

So far we have dealt with the alignment by one component 
of RATs, and pointed out that Qei acts to align grains with 
long axes _L to k, while Qe2 can produce the alignment for 
some grains with long axes || to k. In reality, grains are 
driven simultaneously by both components, so the align- 
ment of grains depends on which component of RATs is 
predominant. 

As discussed in previous sections, for AMO, we have the 
stationary point which is determined by equations (|52|l and 
(|53p . Also, we have analytically shown that the stationary 
point Qa = {n for left helical AMO) is the repellor point. 
To test our expectation, we produce phase trajectory maps 
of grains driven by RATs with the precessing, spinning-up 
and aligning components given by equations (|36|l - (|38|l . Fig. 



° We note that the angle between J and k or B is shown rather 
than the angle between ai and k or B as in DW97. 
^r We checked that even the addition of Qei with an amplitude 
of 10~^ of the amplitude of Qe2 can destroy such an alignment. 



+100 
+ 50 



a 




+ 50 



+ 100 

-1.0 -0.5 0.0 0.5 1.0 

COS© 



AMO, right helicity, B = 0, Q^^ = Q 




Figure 12. Phase trajectory maps of grains in the absence of 
magnetic field corresponding to Qe2 = ( Upper Panel) and 
Qei = {Lower Panel). The map in upper panel exhibits two 
attractor points A, B corresponding to a perfect alignment as 
pointed out by the analysis; E is a low-J attractor point and two 
repellor points C, D. The lower panel shows two repellor points 
A, B and two attractor points C, D in which D corresponds to 
"wrong" alignment. 

[T3]shows the maps for right and left helical AMOs. It can be 
seen that the phase maps have two stationary points which 
positions are exactly given by analytical expressions (|52|l 
and (|53|) . Among these stationary points, there are two re- 
pellor points A, B at Os = (or n for lower panel), and 
one low-J attractor point C at Oj, = tt (or = for lower 
panel). In addition, we see that for the right and left helical 
AMOs, the trajectory maps are related through a transfor- 
mation: O — > TT — O. 

We have seen that , in both cases of grain alignment by 
one component and joint actions of components, the impor- 
tant features (e.g., attractor and repellor points) present in 
the trajectory maps that are constructed with exact RATs 
are consistent with the predictions based on the approximate 
formulae using the fitting function f„/2 and g. It indicates 
that, though the fitting functions do not give the best fit in 
some particular angles 0, they can provide us, intuitively, 
the alignment property of grains by RATs. 

We found in the case of alignment by one component 



AMO, right helicity, B = 0, 0^2 = 
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Figure 13. Upper and lower panels show the phase trajectory 
maps for grains with right and left helicity, respectively. Figs show 
that grains have a perfect alignment with respect to k with one 
low- J attractor point C at = tt (upper panel) and © = (lower 
panel). Also, two repellor points A, B are shown as predicted in 
our analysis. The maps for upper and lower panels are mirror- 
symmetric. 



that if only Qei is at work, it aligns grains with two attractor 
points with J 3> Jth, = 0, i.e., high-J attractor points. 
In contrast, when only Qe2 acts, those points are repellor 
points; Instead Qe2 produces an attractor point at = n/2. 
When both components act simultaneously, are there high-J 
attractor points? Obviously, we may conjecture that if Qei is 
predominant over Qe2, then high-J attractor points should 
still appear. Otherwise, high-J stationary points are repellor 
points. 

For our default AMO with a — 7r/4, we have pre- 
dicted that the stationary points = are always repellor 
points. This can also be understood in terms of the ratio of 



Q 



max / /^max 



1, for this case, i.e, the dominant criteria is 



not satisfied (see equations 1191 and I20|l . 

To study when we have high-J attractor points, let us 
write RATs for AMO in a simplified form 



Qei 
Qc2 



{e2 



-(5 cos^0-2), 
sin 20, 



(61) 
(62) 



where QS^^iQS"^ are maximal values of Qe\,Qe2- 

Clearly, stationary points for this model are 0s = be- 
cause aligning torque -F(0s) = —Qei sin 0s + (5e2 cos 0s = 
at 0s = 0. 

Using the criteria of an attractor point (i.e., equation I55|) . 
it follows that 0s = is an attractor point if 

1 dQe2 , 



-(0s = 0) < 1, 



(63) 



Qei dO 

Plugging equations (|61|l and (|62|l into equation (|63[) . we get 
Q7i" > 2QT2'""'. (64) 



From Fig. IB6l and equation (|64p . it follows that for the orig- 
inal AMO, i.e., AMO in which the relative magnitude of 
RAT component is not rescaled yet, the maximum of the 
ratio maa::[Q™'"/Q™'"] is 1.3, and therefore the stationary 
points 0s = (corresponding to high J) are always repellor 
points. However, as we have discussed earlier, their relative 
magnitude is adjustable. Thus, AMO can produce the phase 
map with high attractor points for ip — (P , provided that 
it satisfies equation ((64}. Note, that the criteria (|64p is only 
applicable for the case of alignment with respect to k or 
■(/' = 0" (i.e., k||B). For an arbitrary angle ip, the criteria is 
shown in the lower panel of Fig. 1241 There it can be seen 
that for some irregular grains (i.e. shape 1, ISRF) 0s = 
can correspond to high attractor points. 

6.7 Alignment for irregular grains 

Similar to the AMO case, we consider first the alignment of 
dust in respect to the radiation direction k. In particular, 
to compare the action of RATs for AMO and an irregular 
grain, we consider the alignment that is induced by individ- 
ual torque components. 

To calculate the phase map, we use again the parame- 
ters from Table [2] We exemplify the alignment for irregular 
grains using shape 1, which shows the maximum deviations 
from AMO predictions. 

Fig. [l4]shows the phase trajectory maps of an irregular 
grain (shape 1) driven by RATs calculated by DDSCAT with 
either Qe2 or Q^i taken to be zero. 

We see that these trajectory maps are similar to those 
constructed by RATs from AMO (see Fig. [12} . Indeed, for 
Qei = case, the map for AMO in Fig. [12] shows two at- 
tractor points at cos = 1, which are found in Fig. 1141 
However, in the situation when only Qe2 acts, the upper 
panel in Fig. 1141 shows an attractor point at = 7r/2 and 
J ^ Jth, which is somewhat different from Fig. 1121 This dif- 
ference stems from the fact that Qe2 for shape 1 is though 
small, but not equal to zero at = 7r/2. Therefore, RATs 
are still able to spin up grains and align them there. 

Fig. [15] shows the phase map when all torque compo- 
nents are at work corresponding to shape 1 and ISRF. It is 
shown that the trajectory map does not have the attractor 
point at = 7r/2. Thus, similarly to AMO, the attractor 
point corresponding to aligning grains with long axes par- 
allel to the direction of radiation disappears when non-zero 
Qei is accounted for. 

Moreover, from Fig. [15] (shape 1) and the upper panel 
in Fig.[T3]for AMO, it is seen that both maps have the same 
repellor point at cos = 1 in the upper panel. However, 
the point C in Fig. [15] is an attractor point, rather than a 
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Figure 14. Phase trajectory maps in tlie absence of magnetic 
field corresponding to Qe2 = ( Upper Panel) and Qei = {Lower 
Panel). The map in the Upper Panel has three attractor points 
A, B, C in which B, C correspond to high angular momentum. 
The map in Lower Panel has one attractor point C corresponding 
to "wrong" alignment and two attractor points with J = 0. 
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Figure 15. Phase trajectory map in the absence of magnetic field 
when all torque components act together shows the alignment 
with one high- J attractor point C , two low- J attractor points A 
and B corresponding to perfect alignment with k, and one repellor 
point D. 




AlisiiimeTH totirdiiniiU; svsleni 



repellor point B as seen for AMO. This difference stems from 
the fact that for AMO, Qei = ■t"'i"i^"2l"2l _L(5cog2Q _ 2) 
is completely symmetric, i.e., Qei(6 = 0) = Qei(6 ~ tt), 
while for the Shape 1, Fig. |6] (the solid-dot line) shows that 
Qei{0 = 0) < Qei{0 = 71"). Therefore, the stationary point 
= in the upper panel, does not satisfy the criteria for 
attractor points, i.e., it is a repellor point, while it is an 
attractor point in the lower panel. 



7 ALIGNMENT WITH RESPECT TO B 

We showed in ij6]that grains can be aligned with respect to 
k. Below we consider the case when magnetic field is essen- 
tial in terms of grain precession (see Fig. I16|) . As earlier, we 
disregard the paramagnetic relaxation. We make an exten- 
sive use of the physical insight obtained with a more simple 
case of alignment in !j6l Indeed, because of the precession 
about magnetic field the analytical treatment of the corre- 
sponding processes gets less transparent here compared to 
that in 



Figure 16. Coordinate systems used to study grain alignment in 
which the external magnetic field B defines the alignment axis. 



T.l Equations of motion in presence of B 

In the presence of magnetic field, equations of motion in 
dimensionless units become 



M 



dt 



dt 



sin ^ 



J 



(65) 
(66) 
(67) 



where f2s is the Larmor precession rate of J around the mag- 
netic field B. Here F,H,G are RAT components projected 
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to three axes ^, j, cf), which are given by (see DW97) 
F{£^,ip,(f)) — Qei(6,0)[— sin tp cos ^ cos (p ~ cos sin £] 

+ Qe2{O,0)[ COS $( COS 1p COS ^ COS (j) 

— sin ?/j sin ^) + sin $ cos ^ sin 0] 

+ Qe3(6,0)[ COS <I> COS 5 sin <;!) + sin "!>( sin i/j sin ^ 

— cos i/j cos ^ cos <;!))], (68) 
H{^,if),(l>) — Qei(6,0)[— sin t/j sin ^ cos 0+ cos t/" cos £] 

+ Qe2(6,0)[ COS $( sin cos ^ 
+ cos ?/j sin ^ cos (p) + sin $ sin ^ sin (^], (69) 
G(5,iA,<?!>) = Q.i(e,0)[ sin i/; sin 4>] 

+ Qe2 [ sin <1> cos (f} — COS "I> cos tp sin <j!>] 

+ Qe3[ COS <1> COS (j) + sin $ COS -i/j sin 0]. (70) 



from to 27r. As a result, 



Here O and <1> are related to ^, ^, <j!> via 



cos — cos ^ cos i/) — sin ^ sin cos 
$ = 



„i sin Q — sin f sin th 
= 2tan ^ : — ^ . \ 



sin ^ sin . 



(71) 
(72) 



Equation (|69|l reveals explicitely that the component 
QesiQ, 0) does not contribute to spinning up grains. On the 
other hand, we found numerically that the last term contain- 
ing Qe3{&,0) in equation (|68|l goes to zero after averaging 
over the precession angle 0. Therefore, similar to the case of 
alignment with respect to k, the only effect of Qe3{Q,0) is 
to induce the grain precession. 

After averaging over the precession angle 0, the equa- 
tions of motion (|65[) - (|67p are reduced to two equations 
for ^ and J, whereas F{^,^, (p), H {^,tp,4)) are replaced by 



7.2 Stationary points for arbitrary shaped grains 

While in this section we deal with AMO, some results can be 
obtained in a general case of arbitrary shaped grains. More 
results of this nature are presented in §8. 

In the presence of magnetic field, aligning and spinning 
torques are complicated functions of RATs, involving ?/;, ^, (j) 
variables. Therefore, it's not easy to derive general analytical 
expressions for stationary points. However, we can find some 
particular physically interesting situations that correspond 
to stationary points. 

For instance, the perfect alignment corresponds to the 
maximal inertia axis ai parallel to the magnetic field, i.e., 
sin = 0. For this angle, from equations (|7ip and (|72p . we 
have = 1/;, and $ = or tt. Hence, Qei(0,O) = Qei{'4',^), 
Qe2{Q,Q) ~ Qe2{il',0)- Equatiou (|68ll becomes 

F{(,s,i',<t>) = Qel(V',0) sin 7/) cos (^+ Qe2(V',0) COS $ COS (j) 

+ Qe3(C,0) COS $ sin (73) 

Obviously, F{^s,'>p, 4>) is a function of the precession an- 
gle (j> about the magnetic field B. Thus, if the grain precesses 
rapidly around B, then we can average F{^a,tp,4') over <j> 



{F{£,s,^))4.^Q.ii^,0) sin / 

Jo 

-FQe2(^,0) COS $ / 

Jo 

■L 



+ Qe3 COS $ 



sin cj>dcj> = 0, 



(74) 



which implies that for grains of an arbitrary shape, and for 
arbitrary direction of light with respect to the magnetic field, 
there are always two stationary points at = and n. 
This very fact makes the alignment of the grains with long 
axes perpendicular to B in some sense the expected one, 
although it does not present a sufficient condition for such 
an alignment. 

If attractor points exist for ^3 different from or n, the 
alignment may get "wrong", i.e. with the maximal inertia 
axis of the grain tending to be parallel to the magnetic field. 
Here and below we adopt the convention that the align- 
ment is "right" if it corresponds to the Davis-Greenstein 
predictions, which made the Davis-Greenstein mechanism 
so popular even in spite of its inefficiency. Needless to say, 
that RATs we seek "right" alignment, i.e. with long grain 
axes perpendicular to magnetic fields without appealing for 
paramagnetic relation. 

The "wrong" alignment may be expected, for instance, 
when the radiation beam is perpendicular to the magnetic 
field, i.e. tp = 7r/2. For this ip, consider the direction of 
"wrong" alignment corresponding to = 7r/2. For tp — it/2, 
and 5 = n/2, we have cos O = — cos (p, and "I> = ■k/2 
for 4> < -K and — 7r/2 for (p > w. Therefore, equation 
becomes 



F(es,V,0) = Qe3(e,0) sin -!>. 



Thus 



(f(C.,^))0= / Oe3(<^,O)d0 = O. 

Jo 



(75) 



(76) 



Here we use the property (5e3(0, 0) ~ sin 2(p in calculating 
the integral. 



7.3 "Right" and "wrong" alignment for AMO 

The introduction of fast precession arising from B makes the 
dynamics of grains more interesting. For instance, it allows 
for a parameter space for "wrong" alignment, high attractor 
points for AMO, shifts of the crossover points. 



7.3.1 Torques considerations 

Above, we found that in the presence of B, two permanent 
stationary points are £,a = 0, n. If being attractor points, 
they correspond to the "right" alignment. 

Consider a case of suspected "wrong" alignment for 
AMO at ~ 7r/2. Since as tp — n/2, cos Q = — cos (p, 
so RATs (see equations 1341 and I35p become 

16Zinin2|n2L(. 2, o\ /'7'7^ 
Qei = ^ — ^(5 COS (p~2), (77) 

40hnin2\n2\ . , n 1 000 2,x ._o\ 

Qe2 = ^-r— ' — - sm 20(1.191 -t- 0.1382 cos (p), (78) 
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Substituting equations (|77p and (|78p into equation (|69p and 
averaging over the precession angle (j>, we get 



16linin2\n2\ 
3A 

40/inin2|n2| 
3A 



(5 cos (p ~ 2) cos 



sin 20 sin <^(1.191)d0 



+ 



sin 2(^(0.1382 cosV)rf<^] = 0. 



(79) 



The fact that the integral (|79p is equal to zero means that 
RATs do not spin up grains when they are perpendicular to 
the magnetic field. 

Now, let us study whether = n/2 satisfies the con- 
dition of an attractor points. Fig. llSl upper shows the spin- 
ning and aligning torques for ip = 89.9" ~ 90" for the 
state ai||J. It shows that, there are stationary points at 

cos ^ = ±1 and cos ~ 0.1, —0.1 for the J parallel and 
anti-parallel to ai corresponding to zeroes of {F{^)}^. For 
the stationary point C with cos = 0.1, Fig. 1181 shows 

^^^^^1^1 > 0, and {H{^s))4, > 0. Hence, this stationary 

point does not satisfy equation (I44p . i.e, it is a repellor point. 
Meanwhile, the stationary point C at cos ~ —0.1, we 



have 



d5 



< and {H{£,s))(f, > 0. As a result, the sta- 
tionary point C' is an attractor point. Similarly, the point A 



< and 



at cos = 1 is an attractor point since — gj- 

{H{(,s))<t> > 0; also, the point B is a low-J attractor point. 
However, the points A' and B' are repellor points in the case 
J anti-parallel to ai. 

Therefore, the possibility of existence of "wrong" align- 
ment is feasible, but it happens at a low angular momentum. 
So, RATs from anisotropic radiation field themselves can not 
maintain the "wrong" alignment with respect to magnetic 
field in the presence of thermal wobbling and the bombard- 
ment by the ambient gas (see Hoang & Lazarian 2007). 

7.3.2 Effect of isotropic torques on "wrong" alignment 

If grains are sufficiently large not to experience frequent 
thermal fiips (see LD99a), they can be subjected to regular 
isotropic torques, which include both Purcell's torques (Pur- 
cell 1979) and those from isotropic radiation fiux (DW96). 
According to LD99b, the this corresponds to grains larger 
than 10~* cm, which is much larger than the typical size 
of the grains in our Table 2. However, such large grains are 
relevant to many astrophysical environments, e.g. comets, 
dark clouds, accretion disks. 

The "wrong" alignment gets modified when we take 
into account RATs induced by isotropic radiation, or equiv- 
alently, Purcell's H2 torques. Indeed, "isotropic" RATs and 
Purcell's spin-wheel are both parallel to the maximal iner- 
tia axis, i.e, ai. Therefore, the total spinning and aligning 
torques are {H{(,,Tp))^ + Qiso, (-F'(C) V'))<?>i respectively. 

Since the aligning torque (F{^,ip))^ depends uniquely 
on RATs induced by anisotropic radiation, the positions of 
"wrong" attractor points do not change. Meanwhile, their 

Here we take 1/) = 89.9" to avoid the singularity that may 
appear due to zero of (H) as shown in equation JTSJ when con- 
sidering condition of attractor points 



angular momentum are added by a term resulted from 
"isotropic" RATs or Purcell pinwheel. 

Now, let introduce H2 torques. H2 torques, after aver- 
aging over the grain rotation around the maximal inertia 
axis, are given by 



I1LOH2 
= ai-;^ p, 



^qas 



(80) 



where 



5 X 10^ ( 4- 



E 



\ 1/2 
0.2eV j 



n{H) 



(81) 



Here is the surface area of a catalytic site of the grain sur- 
face, / is the efficiency of H2 formation, n{H) is the density 
of atomic hydrogen, E is the kinetic energy of the escaping 
H2 molecule, and p is a random variable (see DW97) . In ad- 
dition to parameters given in Table El here we assume that 
n{H)/n = l,f = \,l = 10^°, and E = 0.2eV; we consider 
then the role of H2 torques for p — 1,-1 corresponding to 
Th2 parallel and anti-parallel to ai. 



7.3.3 Trajectory maps 

To illustrate the alignment with respect to the magnetic field 
for AMO, we construct phase trajectory maps using RATs 
obtained by averaging equations (|10p - (|12p . for different ra- 
diation directions ip. To do this we adopt the parameters 
from Table H For 1/; = 0°, i.e., B|lk, the phase map shown 
in the upper panel of Fig. 1 131 exhibits a perfect alignment of 
J with respect to B. 

Fig. [T7]shows the phase trajectory maps for ip = 30, 6O" 
with two repellors A and B, and one low-J attractor point 
C. 

For ~ 90", the lower panel in Fig. [TH] shows an at- 
tractor point C at cos ^ — —0.1 corresponding to J ~ 
and two repellor points at cos ^ = 1, respectively. It is 
shown that the phase map is not symmetric between upper 
and lower panels compared to the maps for ip < 90". This 
stems from the fact that spinning and aligning torques are 
no longer symmetric. The existence of the attractor point 
C indicates that there is, indeed, a "wrong" alignment situ- 
ation as predicted by our analysis above. However, to what 
extend of tp the "wrong" alignment occurs? 

Fig. 1 191 shows the angle of attractor points (upper) and 
their corresponding angular momentum (lower) for a range 
ip = [89.5°, 90°]. In each frame of Fig. Umower. we show 
angular momentum of attractor points present in both upper 
and lower panels of a trajectory map (upper and lower panel 
are labeled). 

For the case Qiso ~ 10~^,p — 0, the upper panel in 
Fig. [19] shows that there are always three attractor points in 
which one attractor point happens at cos ^ = —0.1, —0.05 
and for 1/; = 89.5°, 89.8° and 90°, corresponding to 
"wrong" alignment, two other with cos ^ = ±1. However, 
their angular momentum is very low, about 1.5/iCijt (top 
frame in Fig. \19l ower) . When H2 torques are taken into ac- 
count, if H2 torques are parallel to ai, i.e., p = 1, then the at- 
tractor points cos ^ = ±1 are lifted to J /Iiujt ~ 100, while 
the angular momentum of the "wrong" attractor point is un- 
changed (see the middle frame in Fig. [T9{ouier) . In contrast, 
if p = —1, the angular momentum of the "wrong" attractor 
point can be increased by H2 torques to J / I\lut ~ 100, but 
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100 



AMO, right helicity, ^ = 30^ 



AMO, right helicity, ip ^ 90 




100 



100 



AMO, right helicity, '^=60'' 




100 



Figure 17. Phase trajectory maps for AMO for ip = 30" {Up- 
per Panel) and t/i = 60" {Lower Panel) both have zero attractor 
points C, and two repellor points A, B. 



the attractor points cos = ±1 are unchanged (see lower 
frame in Fig, \19l ower). Note, that the shift of "wrong" at- 
tractor point position toward cos ^ = is also seen on the 
diagram when V ^ 90" fFig. [TOttpperV 

In summary, the presence of "wrong" alignment with 
respect to magnetic fields can be expected, but it happens 
in a very narrow range of in the vicinity of t/j = 7r/2. In ad- 
dition, their angular momentum is very low in the absence of 
isotropic torques such as H2 torques and "isotropic" RATs. 
Naturally, the effects of isotropic torques are negligible if 
grains are flipping, and therefore are thermally trapped as 
it is discussed in LD99a. 



T.4 Alignment for Irregular Grains 

Similarly as in ^we present the phase trajectory maps for 
shape 1, which is the "most irregular" in terms of RATs. 
Consider flrst phase maps for shape 1 obtained for and tp = 
90" first, (see Fig. [20]). For this case, RATs for shape 1 create 
an attractor point at ^ = 90", J /Iilot = 2 in the phase map, 
which corresponds to wrong alignment (see Fig. I20|l . AMO 
also does produce such a "wrong" alignment but at much 
smaller angular velocity, J/Iilut ~ (see Fig. [18}. This 
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Figure 18. Aligning and spinning torques for the particular 
direction of light tp ~ 90" and the corresponding phase map 
for AMO. Upper Panel: Solid line shows the aligning torque 
(i^(5, tp))^ with two zeroes corresponding to two stationary points 
at cos 5 = 0.1,-0.1 corresponding to the case J parallel and 
anti-parallel to ai, besides two zeroes cos ^ = ±1, while the 
dashed line shows the spinning up torque {H(^,tp))^. Lower 
panel: Phase map corresponding to RATs in the upper panel 
shows one "wrong" attractor point C at very low J/Iilut ~ 
and cos ^ = —0.1. 



difference is resulted from the property of Qe2 component. 
Following AMO, Qe2 = at O = 7r/2, while for shape 1, 
this component is not zero at this angle. 

All in all, for various grains shapes studied, we found 
that there is a narrow range -(/> = 85° ^ 90° in which there 
exist "wrong" alignment. 

In addition, we illustrate the alignment in respect to 
magnetic field with phase trajectories obtained for ip = 30". 
For shape 1, similar studies were performed in DW97. How- 
ever, their treatment of crossovers was different from ours. 
Thus the phase trajectories that we observe are different. We 
do not see cyclic maps with grains emerging from a crossover 
to get accelerated by RATs. Instead, we see, similar to the 
case of AMO (see Fig, \17i ivper). the low- J attractor points 
(see Fig[20j. 

Our study provides important insight into the role of 
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Figure 19. Diagram for the position {Upper panel) and angular 
momentum {Lower panel) of attractor points as function of -0 
for three cases, Qiso = — 10~^,p = (no H2 torques), Qiso = 
— 10~^,p = 1 {Tf[2 is parallel to ai), and Qiso = — 10~^,p = —1 
(^H2 is anti-parallel to ai). Each symbol (diamond, triangle, and 
square) denotes an attractor point, and a symbol-line shows the 
extend of the "wrong" alignment as a function of ip- 



Qei and Qe2 components. For instance, in the presence of 
magnetic field Qei tends to provide the alignment with long 
axes perpendicular to magnetic field, i.e. the "right" align- 
ment. The possibility of "wrong" alignment within a lim- 
ited range of (, is related to Qe2- Similarly, the ratio of 
QT\'^ IQT'i'^ determines the existence of the high-J attrac- 
tor points. 



7.5 Ratio of QS'""/Q3"'': existence of high-J and 
shift of low- J attractor points 

The role of the ratio Q™''/Q3'"" on the grain alignment has 
been studied in i]6.6l for the case of no magnetic field (see !j6] 
when this is applicable). Below we consider the role of this 
ratio when the magnetic field is present. We shall show that 
this ratio is important in determining whether grains have 
high-J attractor points. 

In previous subsections we have derived the expressions 
for the stationary points, and presented the trajectory maps 



Figure 20. For i/" = 90", the map shows three high-J attractor 
point A, B and D and two low-J attractor point C and E. Most 
grains in the lower frame of the map align on D with long axes 
parallel to B, i.e., "wrong" alignment. 



100 



ISRF, '(i' = 30°, DDSCAT, shape 1 




100 



Figure 21. Phase trajectory map for grain shape 1 and i^i = 30*^ 
shows the alignment with two repellor points A and B an one low 
attractor point C. 



for the RAT alignment for a particular a, namely, a = 45''. 
It is clearly shown that stationary points cos ^ = ±1 do 
not depend on the magnitude of Qei,Qe2- However, their 
properties, i.e., whether they are attractor or repellor points 
do so. 

Let us first consider a realization of AMO in which the 
functional forms are established for a = 45'' and their mag- 
nitudes are rescaled to have Qri^^/QS"^ = 0.78 (see Fig. 
IB6|) . Note, that this ratio is similar to the ratio of RATs ob- 
tained by DDSCAT for the irregular shape 4 with radiation 
field of wavelength A = 1.2fj,m (see Fig. [5]). 

Fig. [22] shows the obtained trajectory maps for tp — 
60" for both AMO and shape 4. We see that both AMO 
and shape 4 produce the maps with two attractor points 
A, B at cos ^ = —1 and one low-J attractor point C, but 
one difference is that the percentage of grains on A and B 
for shape 4 is higher than for AMO. For AMO and tp — 
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Figure 22. Phase trajectory maps with ?/) = 60 for the case of 
AMO ( Upper Panel) and for the shape 4 (see Fig. |4]l correspond- 
ing to the same ratio Q^^/Qg^^^ = 0.78 {Lower Panel). The 
maps show that both AMO and DDSCAT give rise two high-J 
attractor points A, B at cos ^ = —1 and one low J attractor 
point C. 

60'^, the existence of high-J attractor points A and B with 
QS"7QS"" = 0.78 is not found in the case a = 45° that 
has QS'^^/QS"'^ = 1.2 > 1 in which these points are repellor 
points (see Fig. I17[) . In other words, the existence of high- 
J attractor points depends on the value of QS^'^/QS"'^ as 
predicted in § 16.61 Furthermore, the similarity in the phase 
maps between the AMO and shape 4, which have the similar 
ratio QTi"^ lOT^"^ = 0.78, indicates a good correspondence 
of AMO with irregular grains in terms of grain dynamics. 

While in most cases, we show the phase trajectory maps 
that start from high values of J, in Fig.[23]we also show the 
phase trajectories starting at J = Jth- The dynamics with 
this initial conditions are similar. However, the alignment 
can be achieved faster (see further discussion in i]9.3p . The 
filled arrows in Fig. 1231 mark time scales of O.btgas- We see 
that it takes approximately less than 0.5tgas to get to the 
lower-J attractor point C for grains with initial orientation 
close to C, for the typical interstellar diffuse gas conditions. 
Also, it takes about Stgas to get to the high-J attractor 



AMO, left helicity, f^60° 




100 - 



-1.0 -0.5 0.0 0.5 1.0 

Figure 23. Similar to the upper panel in Fig. 1221 but grains 
start from angular momentum smaller than the value of high-J 
attractor point. Filed circles show initial position of grains. Filled 
arrows show the time interval of O.btgas- Some grains get aligned 
at low attractor point over a short time t < O.btgas, some other 
get to the high J attractor points over t ~ 3tgas- 



points A and B (see Fig. [23]). These time scales are usually 
much smaller than that for paramagnetic damping invoked 
in the Davis-Greenstein mechanism 

Now let us use the approximate functional form of RATs 
given by equations (|61|l and (|62p . and seek the range of 
Q^^^/Q™^, in which the phase map has high J attrac- 
tor points. Fig. shows the ratio of Q;S''^/Q;S"^ for which 
there exist high-J attractor points in the phase trajectory 
map of grains. It is shown that for ip < 45°, Qei is required 
to be dominant over Qe2, at least Q^""^ = to have 

high-J attractor points; their ratio is an increasing function 
of ill. At ip — 45°, one does not have high-J attractor points 
because Qei , Qe2 are equally projected to B and to the di- 
rection ^ perpendicular to B. 

For i/> > 45°, high-J attractor points occur when Qe2 
becomes predominant, i.e., Q™'"/^™'" < 1 (see the lower 
panel in Fig. I24|l . Also in Fig. 1241 the intermediate region 
with parallel lines correspond to the range in which the 
phase trajectory map has only low- J attractor points, while 
high-J stationary points are repellors. 

Fig. [24] shows also the existence of high J attractor 
points for irregular shapes 1 (both monochromatic radiation 
field and ISRF), 2 and 4 (only monochromatic radiation). 
For shape 1 and radiation of A = 1.2fj,m, its phase map has 
low attractor points for < 60° and high J attractor point 
for Ip ^ 60°. But, the map for ISRF has J attractor point 
only in a small range ip < 10°. Interestingly enough, for this 
wavelength, the ratio of exceeds those possible 

for the "original" AMO. As a result, as we mentioned earlier, 
an attractor point is possible for ip — 0, which also describe 
the situation of grain alignment in the absence of magnetic 
field. For such situations AMO has only repellor points (see 

The differences in the range of ip in which the maps 
have high J attractor points predicted by AMO and irreg- 
ular grains (i.e., dot lines extend outside the region sur- 
rounded by thick solid lines in Fig. I24|l exist, but it is mod- 
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f I ) 

Figure 24. Ratio of Q™^/Q™^ for which the phase map has 
high-J attractor points or low-J attractor points as a function of 
ip. Curved solid lines show predictions by AMO, solid lines (maps 
with low-J attractor points only), dot lines (maps with high J 
attractor points) show the result for irregular grains. 
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Figure 25. Upper panel: The shift of position of low-J attractor 
point as function of Q™^/(3™^ for different Lower panel: 
The shift of position of low-J attractor point with respect to -0 
for AMO with three different ratio Q^-^^ IQ^^^ and for irregular 
grains. The shaded area corresponding to the degree of alignment 
R = 1.5 cos^^ — 0.5 ^ shows the range of "wrong" alignment. 



erate. In fact, the correspondence between the ranges of t/" 
for the existence of the high-J attractor point revealed be- 
tween our predictions and the actual irregular grains allows 
one to find out whether the high attractor point is expected 
if only Q^"''/Q™'^ ratio is known. Note, that the existence 
of high attractor points is important for evaluating the de- 
gree of grain alignment. 

Another new effect related to Q^JJ^^'/QeT'' ratio is a reg- 
ular shift of the position of the crossover point, which, as we 
discussed below, is also low-J attractor point. Fig. \2^h vver 
shows this shift as a function of the ratio QS^'^/QS"^ for 
different i/; predicted by AMO. It is shown that the low-J at- 
tractor point tends to shift to cos ^ = as increases. 
Particularly, for ip — 89.9° ~ 90°, the low-J attractor point 
coincides to cos ^ = 0, i.e., grains are aligned with long 
axes perpendicular to magnetic field. The above observed 
tendency is consistent with our earlier discussion in § 16.41 
that Qc2 act to align grains perpendicular to B. 

Fig. \2^lower presents the shift as a function of ^jj for 
particular angles a of AMO and for irregular grains. It is 
clearly seen that, the low-J attractor points for AMO and 
irregular grains are always at cos ^ = 1 for ijj = Qp , which 
means that no any "wrong" alignment possibility exists in 
this case. As ijj increases, they shift to cos ^ = 0, and finally, 
fall on the perpendicular alignment angle cos ^ = at ■(/) ~ 
90". For shape 4, however, the low-J attractor point can 
produce "wrong" alignment earlier than for shape 2, i.e., 
when Ip ~ 85°. This is because the shape 4 has the smaller 
ratio Q™''/Qrf'"' than the shape 2. 

We feel that Q^^'^ jQ^^"^ provides a sufficiently good 
parametrization for torques of irregular grains. Therefore, in 
terms of practical calculations, establishing this ratio may be 
sufficient for describing the alignment of realistic irregular 
shapes. This would require much less computational efforts 
compared with obtaining the shapes for the entire range of 

e. 



8 REGULAR CROSSOVER 

The dynamics of AMO that we observed above was very dif- 
ferent from that of grains in DW97. As we can see in §8, the 
properties of the RAT components are similar for AMO and 
irregular grains, including those studied by DW97. There- 
fore the difference stems from the different treatment of 
crossovers in our and the DW97 models. 

The most striking difference between our and DW97 
trajectory maps is that our crossover points correspond to 
J = and grains cannot get out of these points. Therefore 
our crossover points are also the attractor points. Note, that, 
unlike DW97, we do not observe cyclic maps. The latter are 
the artifact of their model of crossovers adopted there. 

In what follows, we first discuss crossovers in the most 
general terms and later consider a particular case of a 
crossover which is not affected by thermal wobbling (cf. 
WD03) . Then we briefly discuss the possible effects of ther- 
mal wobbling. 

The crossover dynamics has two distinct regimes. The 
first one, described in Lazarian & Draine (1997) takes place 
when the time of the crossover tcros ~ J^/Q where J±_ 
is the angular momentum component perpendicular to ai 
(Spitzer & McGlynn 1979), i.e. the time during which the 
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grain undergoes a regular flipping subjected to the torque Q, 
is shorter than the time of internal relaxation treiax- For H2 
torques Lazarian & Draine (1997) obtained that this regime 
is fulflUed for grains larger that Oc ~ 10~^ cm. When LD99b 
introduced nuclear relaxation it became clear that the typ- 
ical critical size for interstellar grains ~ 10~* cm. This size 
depends on the value of the torques rather weakly, thus we 
may also accept it as the critical size of the grains in diffuse 
ISM in the presence of radiative torques. However, in the 
vicinity of really strong radiation sources the critical size 
gets smaller. Note, that Oc is larger than the "typical" a^ff 
in Table [2j but as we mentioned earlier, grains typical for 
diffuse interstellar gas are different from grains typical for 
other astrophysical environments. 

The randomization of grains during a crossover happens 
due to random processes associated with atomic bombard- 
ment and H2 formation. This is a random walk process in 
which the squared deviation in angular momentum (SJ)^ 
scales with the crossover time tcros- Thus the deviation of 
the angular momentum in the process of a crossover scales 



AMO, right helicity, ^ = 30^ 



5/2 



and 

2 

<=// 



as 5,7/ Jth ~ As Jth ~ {kTahy^^ 

assuming that Q is proportional to the cross-section ~ a 
(see §10), we get the deviation 5J/J ~ '^ijf'^ ■> that decreases 
quickly with the grain size. Therefore the assumption of no 
randomization during crossovers should initially be accurate 
for grains much larger than Oc. Such grains according to Cho 
& Lazarian (2005) are responsible for far infrared polariza- 
tion emanating from dark starless cores. 

For grains less than the critical size Sc LD99a showed 
that the physics is different. As the angular momentum of a 
grain get comparable with Jth, such grain wobbles fast due 
to coupling of rotational and vibrational degrees of freedom 
by the internal relaxation (see Lazarian & Roberge 1997). 
For H2 torques this results in flipping that reverses the di- 
rection of torques and gets the grains "thermally trapped" 
(see LD99a). 

All in all, if grain dynamics is presented in the axes of 
J I Jth and cos 7 with 7 being the angle between ai and J, 
the assumption of J||ai is accurate for J j Jth ^ 1. In the 
range of J j Jth ~ [0, 1] the actual crossover physics should 
be accounted for. 

In view of the discussion above, consider a regular 
crossover first, i.e. assume that a^ff 3> a^. The assump- 
tion that ai is always parallel to J breaks inevitably as 
J|| 0. Indeed, whatever is the efficiency of the inter- 
nal relaxation mechanism, there will be some residual J±. 
Lazarian & Draine (1997) showed that the mean value of 
J± during a regular crossover cannot be smaller than Jth- 
In addition, as J gets small, the efficiency of relaxation drops 
and the gaseous bombardment increases J± further in accor- 
dance with the original Spitzer & McGlynn (1979) theory of 
crossovers. 

To simplify our treatment we disregard all internal re- 
laxation processes during the crossover. This is justifiable 
as for aeff flc the crossover happens on the time scale 
shorter than the internal (e.g. nuclear) relaxation time. Since 
the mirror is assumed weightless, for AMO the dynamics of 
free rotation coincides with that of a spheroid. As a result, 
for a given J, ai precesses around J with a constant angle 7. 
The state of the grain is completely determined by describ- 
ing J in the lab and body systems. Equations of motion for 




AMO, right helicity, ^ = 30'' 




Figure 26. Phase maps for AMO with tfi = 30^. Upper panel 
shows that all grains approach an attractor point C of J = 
corresponding to the same angle cos ^ = — 1. However regular 
crossovers corresponding to Jy = occur at different angles as 
illustrated by lower panel. 



this case are 



Jd cos 7 
di 



£^ = r - 

dt t„ 



dJ 



— cos 7 -I- 

dt 



dJ\\ 
dt 



(82) 



(83) 



where T is RAT, and Jy = J cos 7 is the component of 
angular momentum along the maximal inertia axis. In the 
coordinate system J, ^, 4>, equation (|82p returns to the set 
of equations of motion (|65|) - (|67p . Averaging over the angles, 
"I>, and (j) corresponding to the precession of ai around J 
and J around B, respectively, and solving three resulting 
equations for J, ^, 7, we can obtain trajectory maps during 
the regime of low J and regular crossover. 

In Fig. [21] we show the phase map for the case ip = 30". 
The most striking feature observed there is that the grain 
experiences not a single, but multiple crossovers. Unlike Pur- 
cell's torques, as the grain flips, RATs do not change their 
direction. As a result, they drive J further down. In other 
words, the dynamics of grains roughly corresponds to what 
we observed in the earlier sections of the present paper (cf. 
DW97 where grains are spun up after passing J = Q). 
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Depending on magnetic field value, the crossover may 
happen on the time scale smaller than the Larmor precession 
period. In this case, different grains will undergo crossovers 
at different (j> and therefore will experience different torques 
during crossovers. However, this would not change the qual- 
itative picture, as RATs will drive J to zero irrespectively 
of (p. We shall discuss this issue in more details elsewhere. 

Let us consider qualitatively the issue of thermal fiuctu- 
ations. For the Purcell's torques those can be ignored during 
the crossovers for grains with size a^ff ^ ac- RATs, how- 
ever, are different. They tend to decrease J further rather 
than spin up the grains after a crossover. The time of relax- 
ation, which is also the time over which thermal fluctuations 
of 7 take place, increases as J decreases. Therefore in the ab- 
sence of external collisions one can imagine two situations: 
according to one, RATs stop the grain completely, so no 
thermal wobbling of 7 is possible, another is that J settles 
at the value that is of the order of Jth- The latter seems 
more probable for irregular grains. 

An interesting consequence of the considerations above 
is that a grain tends to get to the J — 0, but the irregularities 
in the motion of the grain axes in respect to the photon flow 
prevents the angular momentum from reaching zero. Thus, 
we can expect that for grains with moments of inertia closer 
to one of a spheroid the value of angular momentum at the 
low-J attractor point is going to be smaller than for more 
irregular grains. 

To find the actual value of effective J, one should recall 
that grains are subject to substantial random excitations in 
the typical astrophysical conditions. Gaseous bombardment, 
interactions with ions, absorption of photons (see Draine 
& Lazarian 1998 for a quantitative description of the pro- 
cesses), affects grain rotation at the low-J attractor point. 
We expect these influences to increase the mean value of J 
above Jth- 

Another effect should also be present in the presence of 
random bombardment. If the phase maps have both high 
and low-J attractor point, the bombardment allow more 
grains to get to a higher-J attractor point. As low-J at- 
tractor points are characterized by higher internal random- 
ization, quite counter-intuitively, the random forcing can in- 
crease the degree of alignment. 

All in all, our considerations above justify the simpli- 
fled model of crossovers that we adopt in the paper. We 
shall provide a more detailed treatment in Hoang & Lazar- 
ian (2007). 



9 PARTICULAR CASES 

Our earlier discussion covered two limiting cases, namely, 
(a) the light beam being the axis of alignment, which cor- 
responds to the precession rate induced by RATs t^^ much 
larger than the Larmor precession rate and (b) when 
the magnetic field constitutes the axes of alignment, i.e. 
tg^ >• f^^. Can this case be ever astrophysically important? 

9.1 Criterion for B alignment 

Consider the precession of ai or J around k driven by the 
component Qes which is perpendicular to the axis ai and k. 



The timescale for RAT precession is defined by 
where 

^ = ^Q,.$ = ^Q.3. (85) 
at Iioj IiLj 

Here cu is the angular velocity of a grain around the max- 
imal inertia axis ai. Qes ~ Qes/lO"^, A = \/l.2fim, 
Urad = Wrad/8.64 X 10~^^ er^s cm~'^, and p = p/Zgcm"^ , 
a_5 = a^f f /I0~^cm, and 7 = 0.1 for anisotropy of ISRF. 
In deriving equation (|84|) the assumption uj — ut is used. It 
is easy to see that for axisymmetric grain shapes, Qei , Qs2 
are equal to zero, while Qes is non-zero (see equations . 
Therefore, the third component produces fast precession of 
grains about k. 

The precession rate above should be compared with the 
Larmor precession rate. A rotating grain acquires a mag- 
netic moment by the Barnet effect which is shown to be 
much stronger than that arising from the rotation of its 
charged body (Dolginov & Mytraphanov 1976). The in- 
teraction of the magnetic moment with the external mag- 
netic field causes the gradual precession of the grain around 
the magnetic field direction. The rate of Larmor precession 
around magnetic field, ts, is given by 

ti3 = 3 X lO'^atsP'^^'^x'^B-'^ s, (86) 

where B = B/S/iG, x = x/3-3 x 10~^ are normalized mag- 
netic field and magnetic susceptibility, respectively. 

The alignment of grains whether with the radiation or 
magnetic field depends on their precession rate around these 
axes. According to equations (|84|) and (|86|) . the ratio of pre- 
cession rate due to radiation and magnetic field is given by 

^ = ^-^^^pxf'-^aZl'i^). (87) 

''B Qe3 Altrad 

It can be easily checked that for typical diffuse ISM, ts/tk ~ 
10"'^, i.e. the Larmor precession is much faster than preces- 
sion induced by radiation, therefore magnetic field is the 
alignment axis. 

For AMO, RATs can be very different as Qe3 arises 
from a spheroidal shape, while the other two components 
arise from a mirror; the ratio of the sizes of the mirror and 
the spheroidal body may be arbitrary. For irregular grains, 
however, we can see all three RAT components to be roughly 
comparable (see Figs I6I8|I . In this case, tk can be used as a 
proxy for the time of "fast alignment" (see i|9.3p . 

9.2 Astrophysical Implications 

It is possible to estimate that for the magnetic field of 5 pG 
grains get aligned with the radiation beam when the density 
of energy in a beam is Uj.g^jj > lO'^ujg-p^p. The appropriate 
radiation fields is typical near stars and supernovae. 

Fig. [27] shows the variation of the ratio of precession 
timescales with distance for different stars in which we use 
Qei = QS""" (for shape 1 and ISRF) for equation (|87)) . It is 
apparent that near the stars, grains precess around k much 
faster than around magnetic field. However, the ratio in- 
creases with distance as the radiation field decreases, so the 
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Figure 27. The ratio of timescales for the precession about the 
hght direction and about the magnetic field assuming that mag- 
netic field is homogeneous with strength B = O.bfiG through 
envelopes. Upper Panel: For three values of star temperature, the 
density dependence on the distance from the star is fixed; Lower 
Panel: For a given star with different functions of density. 
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Figure 28. Alignment in the absence of gas damping for two di- 
rections of radiation ^ = O'' {Upper panel) and 70*^ {Lower panel) . 
Here, the gray circles represent the initial position of grains. In 
the upper and lower panel, the distance between two filled ar- 
rows corresponds to a time interval At = lOtjj, and At = SOtfe, 
respectively. 



precession around magnetic field is faster than around light, 
and grains align with magnetic field at some distance. 

Another case when the alignment can happen with re- 
spect to the direction of light is the case of cometary dust. 
This case may be somewhat more complex, as electric field 
can be present near the comet head. This can induce pre- 
cession of grains with dipole electric moment and therefore 
provide yet another axis of alignment. We discuss this prob- 
lem in Hoang & Lazarian (2007). 

If and tg^ are comparable, the axis of alignment 
does not coincide with either beam or B. This case can be 
relevant to some part of comet grains, but a discussion of it 
is beyond the scope of the present paper. 

Note that is a function of Q. Therefore when its 
amplitude value is larger than tg^ this does not guaran- 
tee that the effect of magnetic field is negligible. For AMO 
the alignment drives grains into the position correspond- 
ing to = it/2, for which the component Qes gets zero 
and therefore the magnetic field again dominates. Thus the 
alignment is expected with respect to an intermediate axis. 



This may be important for explaining circular polarization 
arising from comets (Rosenbush et al. 2007). Naturally, not 
only precession arising from magnetic field, but also from 
electric field and mechanical torques (see §11.4) should be 
taken into account for that case. 



9.3 Fast alignment 

In this section we ignore the damping role of the ambi- 
ent gas, and show that grains can be aligned by radiation 
on a timescale shorter than the gas damping time. Such 
alignment can be called "fast alignment" in analogy with 
the "fast dynamo" process (see Vishniac, Lazarian & Cho 
2003 and ref. therein) that can amplify magnetic field on 
the timescales shorter than the magnetic diffusion time. In 
particular, supernovae flashes can align grains around them 
with respect to the direction of light (see |6]). 

In Fig. l28b pperwe show the map for i/; = 0° in which the 
distance between two arrows represents a time interval equal 
10 precession time tk that is defined in equation ((84}. It can 
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Figure 29. Percentage of grains as function of the ratio of the 
ahgnment timescale due to photon, namely tp^ot to the pre- 
cession time about the radiation direction, tj., for AMO with 
Qmax iQmax _ q (^ppg^ panel) and shape 4 (lower panel). 

Here t = g^am , that is chosen to be equal unity for AMO, and 
e ~ 1 for shape 4. 



be seen that the grains that have initial orientations close 
to the attractor point A require alignment timescales f^hnF^ 
is approximately 40tk to get aligned (see Fig. \28h vver) . In 
contrast, grains that have initial angles far away from the 
attractor point need more time (up to ~ 160tfc) to get there 
(see Fig. [29}. 

Fig. l28l !omer shows a similar effect for the case = 70". 
However, the interval of two arrows represents a time inter- 
val At = SOifc. The map shows that grains bound to the 
low- J attractor point A reach it fast, about 70tk, i.e., on 
the time scales much less than the rotational damping time. 
Some grains that otherwise would go to the high- J attractor 
point stream to infinite J in the absence of damping [lower 
pane/j33- In practical terms the latter fact for this regime 
(i.e., alignment timescale tphot < tgas) is not so important, 



as most grains get aligned at the low- J attractor points any- 
how. We also observe the shift of low- J attractor point A as 
in the upper panel of Fig. 1251 . 

The corresponding time for the fast alignment is pro- 
portional to ratio of grain angular momentum and the com- 
ponent of the torque Trad- In terms of RATs normalized 
components Qr that is related to Trad via equation ((T} 
the relevant combination is given by equation (|84|) . The cor- 
responding function F depends only on Qei and Qe2- The 
amplitudes of those can be very different from the compo- 
nent Qe3 that causes the grain precession. Therefore for such 
grains, the ratio of J /Trad oc J /F can be measured in units 
of "the period of radiation induced precession", namely, in 
tfe oc J/Qes (see equation I84|l . To make AMO more corre- 
spond to irregular grains in this respect, we make the choose 
the amplitude of the AMO components Q^s to be of the sim- 
ilar to the amplitude of Qei- 

Histograms showing the distribution of grains on the 
low-J attractor point as a function of alignment timescale 
corresponding to AMO and an irregular grain shape 4 are 
shown in Fig. 1291 It can be seen that, for AMO, about 
45% grains get aligned with respect to k over i-phot ^ 35ifc 
to 45tfe, and about 22% of grains get there over tphot ~ 
50tfc to 65tfc. A few percent of grains require longer time 
to get aligned, up to 170tfc (see Fig. I29to per') . This relative 
inefficiency of alignment is a consequence of small amplitude 
of the function F in the vicinity of the low-J attractor points 
(F = at the stationary points) . On the other hand, for the 
shape 4, Fig. \29l ower shows that about 55% and 21% of 
grains get aligned with k over 35tk and 65tk, respectively. 
Some others can get aligned over 1704^,. We see that the 
similar distribution of grains as functions of alignment time 
between shape 4 and AMO, though the slight difference in 
percentage of grains corresponding to each tphot /tk present 
due to the fact that, their functional forms of their torques 
are not comletely the same. 

Fast alignment happens in respect to magnetic field pro- 
vided that ts < tk (see equation [ST} , but tphot < tgas- The 



ratio tgas/tB = 1.2 x 10 



a_l which provides a sub- 



stantial parameter space if tphot is much larger that tk, e.g. 
tphot ~ lO^tfc. Comparing tgas in Table 2 and tk given by 
equation (|84[) we may conclude that for typical grains in dif- 
fuse interstellar gas tgas is marginally smaller than tphot and 
therefore the grain phase trajectories are still determined by 

the characteristic time of alignmeniQ- 



gas- However, closer to stars tphot provides the measure of 



10 FITTING FORMULAE FOR RATS 

Astrophysically motivated situations require calculations of 
RATs for grains of different sizes and at many wavelengths. 
This requires rather intensive numerical computations. Our 
encouraging results with AMO motivate us to consider 
whether we can predict the scalings of torques. 



tphot is defined as the alignment timescale of grains due to 
photon or radiation 

We may observe that phase trajectories in Fig. \'2Sl ower di- 
rected to high attractor points also correspond to the aligned 



state of grains, although the stationary state requires tgas to be 
achieved. In this sense all grains get aligned fast. 

We see that the alignment times are more than 30tfc. This 
is due to the fact that the aligning torque gets weaker near the 
low-J attractor points (see Fig. MSh vver) 
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Dolginov & Mytrophanov (1976) associated RATs with 
the scattering of right and left handed photons by a grain. 
For this model one should conjecture that amplitude of 
RATs decreases rapidly with increasing of X/a^ff as the 
grain-photon interactions get into the random walk regime. 
In other words, a sharp peak is expected for the torque 
efflciency for photons with A ~ Oe//- AMO, on the con- 
trary, suggests of a linear increase for A <C le//. Our 
computation, however, indicate that RAT gets constant for 
A < fle//. However, both computational and DDSCAT in- 
trinsic limitations do not allow us to perform calculations 
for A/fle// < 0.1. More studies with other techniques, e.g., 
ray-tracing one, are necessary. 

For grains with A ^ Oe//, Lazarian (1995) suggested 
that the scahng of RATs efficiency should be ~ (A/ae//)~*. 
Such considerations disregard the variations of the optical 
constants. Therefore the testing is essential. 

We calculated RATs in function of X/a^ff for three 
grain shapes. Shape 1 and 2 are shown in Fig. O and a 
hollow grain is produced from the shape 1 by removing the 
core of grain. The latter is done to reduce the amount of 
necessary DDSCAT computations while achieving smaller 
X/a^ff ratios. 

We use both the dielectric function for the smoothed 
astronomical silicate (DW97; Weingartner & Draine 2001; 
Cho & Lazarian 2005) and constant optical constant. Results 
are shown in Figs. I30l and 1311 

We see that the approximate self-similarity (i.e. the de- 
pendence on X/a^ff) is an intrinsic property of radiative 
torques. When optical constant changes as a function of 
wavelength, RAT efficiencies for different grain sizes mostly 
differ at wavelengths corresponding to resonance absorption 
features. 

In addition, RATs have nearly constant magnitude as 
A ~ fle/ / , and decrease steeply with the ratio of wavelength 
to grain size. This is because the scattering of photon by 
irregular grains are strongest as A ~ o-eff- We can fit our 
calculations for RATs (see Figs 1301 13ip by a simple function 
given by 

Qr = 0.4 for — < 1.8, 

= 0.4( — )-" for - > 1.8, (88) 

where a is the spectral index that according to Figs 1301 and 
1311 is between 3 and 4. The first case provides a good fit in 
the whole range of A/oe//. In contrast, the latter case gives 
better fit for the range of A/Oe// < 20. Cho & Lazarian 
(2006) use the former fit formulae to calculate polarization 
degree for accretion disks because grains there are widely 
believed to be very large. 

For the case when ai makes an angle O = 45" with 
respect to k, we found that the self similarity is also valid. 
However, the curve of RATs is shallower, and can be fitted 
by a power index rj = —2. 

To study the efficiency of the self similarity, we calculate 
rotation velocity of grains induced by RATs in which RATs 
are directly computed from DDSCAT and derived from the 
self similarity assumption. We use radiation intensity of a 
molecular cloud (see Mathis 1983) to calculate RATs for dif- 
ferent optical depths A^. Resulting rotation speed obtained 
with these two methods (see Fig. I32p shows clearly that the 
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Figure 30. Upper panel: RATs for grain shape 1, as the function 
of wavelength to grain size for the direction = 0^. We see that 
RATs decrease steeply as X/a^-ff increases. Around X/a^^ff ^ 1 
the dependence flattens. Several positions at which RATs de- 
crease dramatically correspond to the peculiarities of dielectric 
susceptibility. If those peculiarities are ignored, RATs are self- 
similar. Lower panel: The same as the upper panel but for shape 
2. RATs can be fltted by analytical functions with r; = 3 or 4 (see 
equation I88I I. 



self similarity provides a fair agreement between exact cal- 
culations and those based on the self-similarity arguments. 
This allows to reduce the DDSCAT computational efforts 
substantially. 



11 DISCUSSION 

Our study above have approached an important problem 
of the RAT alignment mechanism by studying the funda- 
mental properties of RATs. The goal of such studies was to 
change the status of the RAT alignment from an empirical 
fact to a theoretically understood process. Our work indi- 
cates that grain helicity is an essential property of realistic 
grains, which suggests that it should be accounted not only 
for RATs, but for mechanical alignment as well. 
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Figure 31. RATs for shape 1 and hollow grain with constant 
refractive index: self-similarity exhibits extremely well. They can 
be fitted by a step function and a power function with 77 = 3, 4 
( Upper Panel) and t] = 4,5 {Lower Panel) where t] is defined in 
equation II88I I. 



11.1 Evolution of ideas on RATs 

For the first time grain spin-up due to diflerential scattering 
of left and right handed photons was considered in Dolgi- 
nov (1972). The suggestion was limited to quartz grains, 
however. In Dolginov & Mytrophanov (1976) it was realized 
that grains of made of more accepted astrophysical materials 
can be spun up, provided that they are somewhat twisted. 
This ground-breaking study qualitatively considered both 
the possibility of the RAT alignment with respect to mag- 
netic field and the radiation beam. Regretfully, it was not 
appreciated in its time. 

A reason of the lukewarm response of the community 
to this work was probably related to both the absence of 
reliable estimates of RATs efficiencies and to the fact that 
some features of grain dynamics were unknown at that time. 
For instance, Dolginov & Mytrophanov (1976) claimed that 
prolate grains should align along magnetic field, while oblate 
grains should align perpendicular to magnetic field. The for- 
mer conclusion was based on the assumption that grain rota- 
tion along the axis of minimal inertia is stable. It is only later 
that Purcell (1979) discovered so-called Barnett relaxation 
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Figure 32. Ratio of rotation speed to thermal rotation speed 
with respect to grain size obtained with exact RATs from 
DDSCAT {Upper panel) and using the self-similarity {Lower 
panel). Here, radiation direction is assumed to be parallel to the 
maximal inertia axis ai, that is constrained to be frozen to an- 
gular momentum J. The figures show the rapid increase of the 
rotational speed with grain size, and even in the cloud core with 
high visual extinction Ay = 10, large grain still have suprather- 
mal rotation rate, ensuring that grain alignment is not susceptible 
to random collision by atomic gas. 



and concluded that only rotation about the axis of maxi- 
mal inertia is stableF^. On the basis of this Lazarian (1995) 
claimed that the alignment of both prolate and oblate helical 
grains should happen with long axes perpendicular to mag- 
netic field. However, it was not quite clear at that moment 
what makes the grain helical. In addition, Lazarian (1995) 
underestimated the importance of the RAT alignment. 

The quantitative stage of RAT studies was initiated by 
Draine & Weingartner (DW96 and DW97). There RATs 
were shown to be generic for irregular (and therefore real- 
istic) grains and the magnitude of RATs was reliably quan- 
tified. The alignment with respect to magnetic field was 



The Barnett magnetization that induced Ed Purcell to think 
about the relaxation was described in the same Dolginov &; 
Mytrophanov (1976) paper. 
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demonstrated in numerical simulations. On the basis of these 
simulations DW97 claimed that the RAT alignment was the 
dominant process responsible for polarization arising from 
dust in the diffuse ISM. However, the RAT alignment was an 
empirical fact in DW97 that was based on a limited sample. 
In this situation additional important complexities arising 
from including thermal wobbling of grains (WD03) , gaseous 
bombardment etc. could not help to improving our physical 
understanding of the RAT alignment. 

What we have done in the present paper is we at- 
tempted to clarify the basics of the RAT alignment by seek- 
ing the generic properties of RATs that induce the align- 
ment. We used analytical modeling which was tested with 
numerical DDSCAT simulations. As we discuss further, we 
hope that this work contributes to both more intuitive un- 
derstanding of the alignment and to further elaborating of 
the mechanism in order to get precise predictions of the 
alignment degree in different astrophysical situations. 

It is worth mentioning that even we when know the an- 
alytical form of the torques, the dynamics of the system does 
not get completely trivial. It definitely exhibits interesting 
properties. 

11.2 Our approach 

While RATs were initially treated as a quantum effect aris- 
ing from the difference of scattering of left and right handed 
photons (see Dolginov & Mytrophanov 1976), above we pre- 
sented an entirely classical model of a grain (see Fig. [2| that 
reproduces well their properties (see Fig [5}. 

The gist of our approach above is to consider the basic 
generic properties of RATs and to relate of these properties 
with the RAT alignment. AMO plays a central role in our 
considerations. Our toy model of a helical grain allowed an 
analytical description, which enabled us to treat RATs an- 
alytically. In our study we concentrated on the properties 
of RAT components in the scattering system (see Fig. [T|, 
i.e. on the properties of Qei, Qe2 and QeS. These compo- 
nents show a remarkable similarity for grains of very dif- 
ferent shapes (see Fig. and AMO (see Fig. [S}. The chi 
squared test we present in i^4.5l returns the mean value of 
for both Qei and Qe2 about 0.2 (see Fig. I10|). This provided 
us with the empirical justification of using AMO for studies 
of RAT alignment. AMO provides us with both useful intu- 
itive model to think about the alignment and analytical for- 
mulae that allow straightforward quantitative calculations. 

We find that the basic properties of RATs obtained with 
AMO are very similar to the basic properties of RATs for 
arbitrary shaped grains. This, for instance, allows us to talk 
about helicity of grains being the most important attribute 
for the RAT alignment. 

Our major goal above was to get a better understand- 
ing of the physics of the RAT alignment. To do this we 
adopt a model, similarly to one in DW97, that disregards 
the wobbling of the grain axes with respect to angular mo- 
mentum direction (cf. Lazarian 1994; Lazarian & Roberge 
1997; LD99ab), but treats crossovers differently, i.e. in the 
spirit of the Spitzer & McGlynn (1979) model. This provides 
a substantial change in the dynamics of grains. For instance, 
we do not observe cyclic maps reported in DW97. 

While studying the properties of RATs we addressed 
the question of the necessary conditions for RAT alignment 



to happen and to fail. We also study the RAT alignment 
that takes place in the radiative-dominated environments, 
where the direction of radiation defines the axis of alignment 
(see Dolginov & Mytrophanov 1979; Lazarian 2003). Such 
an alignment is both astrophysically important and provides 
a good insight into the physics of the RAT alignment. 

11.3 Accomplishments and Limitations of the 
present study 

We feel that our major accomplishment in the paper above 
was establishment of the analytical form of RATs and clar- 
ification of the role of different RAT components. We hope 
that AMO unveils the mystery that have surrounded the 
RAT alignment from the time of the mechanism introduc- 
tion. 

Another important conclusion that follows from AMO 
is that the RAT alignment is not limited to grains with A ~ 
a^ff, as it was believed before. As the similarities between 
the torques that we obtained for AMO in the A ^ a^f / limit 
and for irregular grains in the A > fle// limit are striking, 
our work shows that the RAT alignment should take place 
also for large grains, that are present in accretion disks and 
dark cloud cores (see Cho & Lazarian 2005, 2007). 

Obtaining generic properties of RATs makes the RAT 
alignment more a predictable theory and opens avenues for 
further theoretical advances, e.g. including thermal fluctua- 
tions, random bombardment, H2 torques etc. Our establish- 
ing of a subdominant nature of one of the RAT components, 
namely Qes, simplifies the theoretical treatment of RATs. 
Insights into the generic properties of two other components 
allow to reduce the amount of numerical computations nec- 
essary to determine the degree of achievable alignment. For 
instance, in the current paper, we found that the existence 
of high- J attractor points depends on the Q^^'/Q™^ ra- 
tio. For practical applications it is important, that the cri- 
terion for this established with AMO works well with ir- 
regular grains. We note parathentically, that while Qei and 
Qe2 demonstrate universal behavior, the aligning and spin- 
ning torques (see equations (|68p and (|69p ) that present their 
combinations do not demonstrate this. 

The present study identifies a parameter space for which 
the RAT alignment may be suspected to be "wrong", i.e. to 
happen with long axes parallel to magnetic field. In addition, 
it provides simple scalings for RAT dependences on the ratio 
of the radiation wavelength to the grain size. 

Our study reveals new properties of the RAT alignment. 
First of all, the alignment may be fast, i.e. happen in a small 
fraction of gas damping time (see ij9.3p . This has important 
consequences for the environments with fast changing radi- 
ation, i.e. circumstellar regions, interstellar medium in the 
vicinity of supernovae fiashes etc. Moreover, we could see, 
that the alignment is different depending whether the initial 
angular momentum is small or large (see Figs 1221 and I23p . 

The approximate self-similarity of RATs (see §10) is an- 
other practically useful property of RATs. Combined with 
the established universality of the functional form of the 
components Qei and Qe2 and the established dependences 
of the properties of trajectory maps on the ratio of these 
components, this allows radically reduce DDSCAT calcula- 
tions that may be necessary to find the expected degree of 
alignment for an ensemble of realistic grains subjected to a 
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realistic radiation field. In fact, we find that RATs change 
the ratio of Q^^'/QS"^ with A and this is the most impor- 
tant difference that the variation of the radiation wavelength 
entails. The functional forms of the torques do not change 
much and can be well approximated with those of AMO. 

One limitation of AMO is the upper limit of the ratio 
QS^'^/QS'"', that makes AMO more appropriate to irregu- 
lar grains with A > 3ae// and A < Oe//. It indicates that, 
though AMO is established based on the geometric optics, 
i.e., A ^ cieff, it is also applicable for the opposite limit. 

In more general terms, our study proved that irregular 
grains can be characterized by helicity. Grain rotation pro- 
vides the averaging that defines the helicity axis, while the 
irregularities define whether the helicity is left or right. The 
phase trajectories of grains that are the mirror-symmetric 
images of each other are mirror-symmetric (see Fig. I13p . As 
expected, the torque component Qes, unlike the other two 
components, coincide for an irregular grain and its mirror- 
symmetric image (see the lower panel in Fig. [S}. Indeed, 
this component is subdominant for most of the alignment 
processes and not related to grain helicity. 

In our study we do not directly address the grain align- 
ment efficiency. Some statements can be made, however. For 
instance, when grains are aligned rotating suprathermally 
the direction of J is immune to the randomization arising 
from the gaseous bombardment. In addition in this case, 
J||ai provides a good approximation. We find, however, that 
an appreciable subset of grains rotates with thermal veloc- 
ities. For those the internal randomization of grain axes in 
respect to J may be important. Does this signify a new crisis 
of the grain alignment theory? We do not believe so. Even 
in the absence of high- J attractor points RATs will drive J 
back to low attractor points, which in most cases, as we dis- 
cussed in the paper, correspond to the preferential alignment 
of grains with long axes perpendicular to magnetic field. As 
for the internal alignment, according to Roberge & Lazarian 
(1997) for typical interstellar grains, this alignment is tan- 
gible even for J ~ Jth- A detailed study of the attainable 
degrees of alignment is provided elsewhere. 

We have not discussed RATs of the strongly absorbing 
materials, e.g. graphite. We expect the torque components 
to show more irregularity for such grains. As the grain align- 
ment theory matures and extends to the environments dif- 
ferent from molecular clouds and diffuse interstellar gas (see 
Lazarian 2007), the importance of the studies of wider range 
of materials will get more pressing. 

11.4 Rates of Alignment and Rotation 

As we mentioned above the RAT alignment can happen on 
time scales much shorter than the gaseous damping time. 
This finding corresponds to the notion in Dolginov & Mytro- 
phanov (1976), that the alignment happens on the time scale 
that is required for the radiative torques to deposit a grain 
with the angular momentum of the order of its initial an- 
gular momentum. Such a fast alignment makes grains good 
tracers of magnetic field when radiation direction changes 
quickly. 

The fast alignment takes place for low- J attractor 
points of the grain phase trajectory map. These are the most 
probable attractor for the grains to end up with. Thus, most 
grains do not rotate suprathermally when subject to RATs. 



In this sense the RAT alignment tends to minimize grain 
angular momentum. 

While we expect that in the presence of thermal wob- 
bling and gaseous bombardment most grains will rotate 
thermally, there is a radical difference between this effect 
and the effect of thermal trapping discussed in LD99ab. The 
effect of thermal trapping there is based on the compensa- 
tion of the Purcell rocket torques, e.g. those related to H2 
formation, by thermal flipping of grains. The more efficient 
thermal flipping, the more efficient is the trapping and the 
less chance of a grain to get high angular momentum. On 
the contrary, we have seen in ^ and [7| that without thermal 
wobbling, the significant fraction of grains ends up in the 
state of J = 0. In other words, thermal fiuctuations increase 
the value of J to a higher, i.e. thermal value. 

In spite of the fact, that most of the grains tend to ro- 
tate with velocities much less than the maximal velocities, 
iOmax, that RATs can spin the grain up, we believe that 
the parametrization of the alignment in terms of uimax/u^T, 
where ujt is the thermal rotational velocity, may be a rough 
practical way of describing alignment. Indeed, the above ra- 
tio reflects the relative importance of RATs compared with 
those related to gas. When RATs force the grain into a low- 
J attractor point, their ability to do this would also depend 
on this ratio. A further research should reveal more sophis- 
ticated and precise parametrization of the RAT alignment, 
however. This parametrization is necessary, for instance, to 
predict the expected alignment from the numerical simula- 
tions of magnetized molecular clouds (see Cho & Lazarian 
(2005); Pelkonen et al. (2007); Bethell et al. (2007)). 

We have discussed in H7.5I that for some phase trajec- 
tories high attractor points are available. The suprather- 
mally rotating grains correspond to high- J attractor points. 
It takes them about 3 damping times to reach such points 
(see also DW97). However, our analysis shows that grains 
get aligned even before they reach high-J attractor points. 
Therefore the RAT alignment can happen over shorter time 
scales for all grains provided that the radiation is intensive 
enough. 

The predominance of low- J attractor points has conse- 
quences that go beyond the problems of grain alignment. If 
grains rotate slowly, then loosely connected conglomerates 
constituting fractal grains can exist. Ever since the classical 
work by Purcell (1979), the suprathermal rotation had been 
thought to destroy such grains. When LD99ab showed that 
Purcell's torques may not be capable to spin-up grains less 
than 10~* cm, it was still thought that radiative torques can 
do the job. Our work questions this (see also WD03) . 

11.5 Direction of Alignment 

The alignment may happen with respect to radiation rather 
than to magnetic held if the precession induced by radiative 
torques is faster than the Larmor one (see !j9|. We found, 
that in the presence of magnetic held the alignment can still 
happen in respect to the direction of the beam or, equiva- 
lently, the direction of the anisotropy of radiation, provided 
that the rate of precession arising from the radiative torques 
is faster than rate of the Larmor precession. This is the case 
of comets sufficiently close to the Sun, ISM in the vicinity 
of supernovae and some circumstellar regions. Over vast ex- 
panses of diffuse interstellar medium and molecular clouds. 
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however, the generic RAT alignment is with respect to mag- 
netic field, thus enabling easy tracing of magnetic fields via 
polarimetry. Note, that our study shows that many features 
characteristic of the alignment in the absence of magnetic 
field carry over to the case when magnetic field is present. 

Our important finding is that while the generic align- 
ment is "right", i.e. with the long grain axes perpendicular 
to magnetic field, for a range of angles ip between the mag- 
netic field and the direction of the beam around = n/2, 
the alignment may be "wrong", i.e. it happens with long 
grain axes parallel to magnetic field. However, the range of 
the angles is rather narrow. As a result, we do not expect 
the effect of "wrong alignment" to persist when grains un- 
dergo thermal wobbling. This wobbling is likely to vary the 
direction of the grain axes with respect to the light direction 
beyond the range angles in which the alignment is "wrong" . 



11.6 Magnetic Field and Gas Streaming 

Unlike DW97, in the paper above we disregarded the ef- 
fects of paramagnetic alignment altogether. When we con- 
sider dynamically important field, its only effect is to induce 
averaging due to Larmor precession. We believe that our ap- 
proach is correct, as for paramagnetic grains the effects of 
paramagnetic relaxation are marginal on times over which 
the RAT alignment takes place. 

Gas streaming can induce its own alignment direction. 
Dolginov & Mytrophanov (1976) assumed that magnetic 
field or a gaseous fiow defines the axis of alignment depend- 
ing on the ratio of Larmor precession time to that of me- 
chanical alignment. On the basis of our study of Qe3 with 
AMO, we believe that a more physically motivated distinc- 
tion is related to the rate of Larmor precession versus the 
precession arising from the mechanical analogy of the Qe3 
torque. Such torque for an spheroidal grain can be obtained 
from formulae in Appendix B by substituting the value of 
the gas atom momentum mv instead of the photon momen- 
tum. The corresponding precession timescales ratio is (see 
equation I85|l 

, S,"Ti-0-5 ^-0.5 £, 

^=3x10^ ^^^, "-^ . , (89) 

where Qe3,gas is the third component of torques induced 
by the gaseous flow, which is the analog of Qes for RATs. 
In equation (|89p . a — ci x 0.1 is the probability of elas- 
tic collision, Vfion, = Vfiow X vthermai is gas flow vclocity, 
nn ~ fiH X 30, T = T X 100 is gas density and temperature, 
respectively, equation (|89|) indicates that for sufficiently in- 
tensive gaseous flows, for instance, Vfiow > 10^, the align- 
ment will indeed happen with respect to the flow direction. 
Note, that we predict that gaseous flows may define the di- 
rection of alignment for a wider parameter range compared 
to that in Dolginov & Mytrophanov (1976). Moreover, we 
claim that mechanical flows can define the axis of alignment 
even for subsonic flow velocities, i.e. at those velocities for 
which the Gold alignment and its modifications (cf. below, 
however) are marginal. 

A conceivable situation is that the gaseous bombard- 
ment arising from grain streaming defines the axis of align- 
ment, while RATs do the alignment job. This situation takes 



place when Qe3,RATs is less than Qei,gas, which defines 

tflow "^rad 10^ ^rad f90) 

tk ^ OLmHv)i^^nH ~ av)i^^fhH ' 

where Urad = Wrad X UISRF- 

This can be the case of alignment in a part of comet at- 
mospherj^. Naturally, combining equations (|89p and (|90|l it 
is possible to establish when streaming defines the alignment 
axis in spite of the presence of magnetic field. 



11.7 AMO and Mechanical Alignment of Helical 
Grains 

Our present paper is devoted to RATs and the alignment 
that they entail. However, our consideration of a helical 
grain is quite general. In fact, the functional dependence 
of the torques that we obtain for our model grain is valid 
when atoms rather than photons are reflected from the mir- 
ror. Therefore we may predict that for elastic gas-grain col- 
lisions the helical grain^^ will align with long grain axes 
perpendicular to the fiow in the absence of magnetic field 
and with long axes perpendicular to B, when dynamically 
important magnetic field is present. If atoms attach to the 
grain surface and then are thermally ejected from it, this 
changes the values of torques by a factor of order unity. The 
only way that the uncompensated torques can vanish for a 
helical grain is if the correlation is lost between the place 
at which the atom strikes the grain surface and leaves ilF^. 
The latter is rather improbable for sufficiently large grains. 
As the result, one has to conclude that helicity is a generic 
property of the interaction of irregular grains and atomic 
fiows. 

This conclusion alters substantially the paradigm of 
grain motions in diffuse gas. Since the time when Gold 
(1951) proposed his first simple model of mechanical align- 
ment, it was considered essential to have supersonic grain- 
gas velocities to achieve any tangible alignment. Indeed, 
both the Gold original process and the ones proposed for 
suprathermally rotating grains, namely cross-sectional and 
crossover alignments (Lazarian 1995) require the supersonic 
drift to ensure that the momentum deposited in a regular 
way exceeds one deposited due to thermal atomic motions. 
This is not a requirement for the alignment of helical grains! 
For those the regular momentum is deposited in proportion 
to the number of collisions, while the randomization adds 
up only as a random walk. In fact, the difference between 
the mechanical alignment of spheroidal and helical grains 
is similar to the difference between the stochastic Harwit 
(1970) alignment by stochastic absorption of photons and 
the RAT alignment. While the Harwit alignment requires 



Another case also relevant to the comet atmosphere is that 
the grains have electric dipole moments, while comet atmosphere 
has electric field. Then the comet electric field defines the axis of 
alignment. 

The mechanical alignment of helical grains was briefly dis- 
cussed in Lazarian (1995) and Lazarian, Goodman & Myers 
(1997), but was not elaborated there. 

Even in this case the local anisotropies of the surface at the 
place of atom impact can result in efi'ective helicity similar to the 
case of damped oscillator in Fig. 1181 (lower panel) 
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very special conditions dominate, the regular RATs easily 
beat randomization. 

Similarly as in the case of RATs, where it is frequently 
possible to disregard the Harwit process, it should be possi- 
ble to disregard the Gold-type processe^^ for the alignment 
of helical grains. Note, that because of the property of he- 
licity not to change sign during grain flipping, we do not 
expect to observe the thermal trapping effects described in 
Lazarian & Draine (1999a) to be present for the mechani- 
cal spin-up of helical grains. The effects that decrease the 
efficiency of the mechanical alignment of helical grains are 
discussed elsewhere. 

In typical conditions of diffuse interstellar medium, the 
mechanical alignment of helical grains tends to act to align 
grains in the same direction as the RATs, i.e. with longer 
axes perpendicular to magnetic field. The relative role of 
the two mechanisms should be revealed by further research. 
The currently available data (see Lazarian 2007) agrees with 
the RAT mechanism being the primary source of alignment. 
However, the situations are possible, when mechanical align- 
ment reveals magnetic field, when RATs fail to do so. 



12 SUMMARY 

In this paper, we studied the properties of RATs, and 
how different RAT components affect the grain alignment. 
Briefly, our results are as follows: 

1. We found that a simple model of a helical grain which 
consists of a reflecting spheroidal grain with an attached 
mirror reproduces well the functional dependences of RATs 
obtained for irregular grains using DDSCAT. 

2. From the generic properties of RATs we predicted the 
preferential alignment of grains with long axes perpendicular 
to the direction towards the source of light, provided that 
magnetic field effect is subdominant. 

3. The magnetic field is important and defines the axis 
of alignment when it induces the Larmor precession that is 
faster than the precession arising from the Qes-component of 
RATs, i.e. the component not related to grain helicity. This 
component is present for spheroidal grain, for instance. 

4. When magnetic field is important, RATs tend to pro- 
vide both the alignment of long grain axes perpendicular 
to magnetic field. With or without magnetic field, most of 
grains are driven to the low angular momentum attractor 
points, as RATs align grains. Grains can be driven to the 
low- J attractor points on the time scales much less than the 
gaseous damping time. The very existence of the high angu- 
lar momentum attractor points, and therefore grains rotat- 
ing much faster than the thermal velocity, is not a default 
and depends on the ratio Q^i and Qe2 components. 

5. RATs can also induce "wrong alignment", i.e., the 
alignment with long grain axes parallel to magnetic field. 
The range of angles for "wrong alignment" is narrowly cen- 
tered around the 7r/2 angle between the direction of light 
and magnetic field. This range is expected to be diminished 

This alignment tends to minimize grain cross section, which 
means, for instance, that for grains streaming along magnetic 
fields the non-helical stochastic torques will tend to align grains 
with longer axes parallel to magnetic field, while helical torques 
will tend to align in the perpendicular direction. 



when thermal fluctuations are accounted for. Thus the RAT 
alignment is capable to account for most of the observed 
polarization. 

6. RATs exhibit approximate self-similarity that allows 
one to be expressed them as a function of the ratio of the 
wavelength to grain size. The dynamics of the grains can be 
reproduced relatively accurately when the self-similarity is 
used. 

7. The RAT alignment is a particular case of the align- 
ment of the helical grains. Therefore our results can be gen- 
eralized to describe the mechanical alignment of irregular 
grains. Such an alignment is efficient for both supersonic 
and subsonic gaseous flows. The mechanical alignment may 
happen either with respect to the magnetic fleld or the di- 
rection of the flow depending on the rates of precession that 
are induce by the flow and magnetic field, respectively. 
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APPENDIX A: REFLECTING OBLATE 
SPHEROID 

The toy model of a grain (see Fig. [2} allows us to derive 
analytical formulae for RATs. Let us consider RATs for 
an oblate reflecting spheroid characterized by semi-axes a, b 
where s = a/b < 1 (see Fig. IA1|) . ai,a2,a3 are principal 
axes of the spheroid with moments of inertia Ii > I2 ~ I3, 
respectively. Assuming that a photon beam of wavelength A 
is shined along k||ei, the photon reflections happen on the 
spheroidal surface at the location determined by a normal 
unit vector N and a radius r. For the sake of simplicity, a 
perfect reflection is always assumed in this study. The loca- 
tion of impact on the grain surface is specifled by the radius 
r and the normal vector N, which are, respectively, given by 

r — a sin ryai + b cos rj cos ^a2 -I- 6 cos 77 sin ^£13, (Al) 

N = ai sin rjai + 61 cos i] cos ^a2 + &i cos r; sin ^£13, 

(A2) 

where ai = [ sin ^j?-|-(l-e^ cos ^77]"^''^, 61 = ai(l-e^)^/^, 
and e is the eccentricity of the spheroid, ^ — [0,-K],rj = 
[-7r/2,7r/2] (see Fig.|Xl] see also Roberge et al. 1993). 

Due to the symmetry around ai, we only need to find 
RATs for a single rotation angle, e.g. for (3 = 0. Therefore, 
when the grain axis ai makes an angle O with respect to 
the photon beam, we have 

ai = cos Oei + sin Oe2, (A3) 
a2 = — sin Bei + cos Oe2, (A4) 
a3 =63. (A5) 

Substituting equations HA3|) - (|A5|) into (|AH) and (|A2|) we ob- 
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Figure Al. Upper panel represents the coordinates and vectors 
for tlie spheroid. Lower panel shows the fitting function K{0, e) 
for Qe3, depending on the eccentricity of the spheroid e and its 
angle with light beam ©. 



tain 

r — {a cos O sin rj — b sin O cos rj cos ^)ei 
+ (a sin O sin + 6 cos Q cos rj cos 5)e2 
+ b cos rj sin ^63, (A6) 
N = (oi cos O sin r) — bi sin cos rj cos ^)ei 
+ (oi sin Q sin rj + b\ cos O cos rj cos ^)e2 
+ &i cos 77 sin ^63. (A7) 

Hence, RAT produced by the photon beam is defined by 

dTrad = 7r X AP = -2p7FdA(k.N)[r x N], (A8) 

where 7 is the anisotropy degree of radiation field, pph is 
the momentum of each photon, F is the flux of the incident 
hght beam, and dA is an area element on the grain surface 
given by 

dA = eb^f{rj) cos rjdrjdi, (A9) 



1— p2 , . p 



Substituting 



hc/\ 



and dA 

h ' 



from equation (|A9P into the above equation, we obtain 

7MradA6^ 4e 



-(ei.N)[(r2A''3 -r3Af2)ei 



2 ^ A 

+ (r3iVi - riA''3)e2 + (ri7V2 - r2A''i)e3]/(j7) cos rjdrjd^. 

(All) 

From equations HA7|) and (|A8|) . we get 



r-2A'"3 - r3iV2 

r-37Vi - riA'^s 
riN2 — r2Afi 



abi 




bai 




2 




abi 




bai 




2 




abi 




bai 



sin 2r? sin ^ sin 6, (A12) 
sin 277 sin ^ cos Q, (A13) 
sin 277 cos 5. (A14) 



Plugging equations (IA12|) - (IA14|) into equation (|Alip . and 
integrating over the surface, we obtain 



i rad = ^ (yeiei + Qe2e2 + QeSGS) 

where the RAT components are given by 

_4e abi - bai 



(A15) 



Q 



A 



\ai sin 77 cos 



7r/2 ^0 

61 COS 77 cos ^ sin B| sin 277 sin ^ sin 0/(77) cos rjd^drj, 

(A16) 



A 



|ai sin ?7 cos 



tt/2 Jo 

bi cos 77 cos ^ sin Oj sin 277 sin ^ cos 0/(77) cos rjd^drj, 

(A17) 



ir/2 JO 



A 



|ai sin ?7 cos 



— 61 cos 77 cos ^ sin Oj sin 277 cos £,f{rj) cos rjd^drj. 

(A18) 

Note, that we use the absolute value for ei.N because in the 
integral over 77, we always use the range [—n/2,n/2]. 

Due to the presence of the term sin ^ in equations 
(|A16|) and (|A17p . their integrals over the range ^ — [0,n] 
vanish. The integral given by equation (IA18|) provides us a 
function of O which can be fitted by a function sin 20 and 
a fitting factor 



?e3(e) = ^(s^ - l)K{e,e) sin 26, 



(A19) 



where /(r7) — 
Therefore 

dVrad = -2pphjF{ei.N)dA[{r2N3 - r-3iV2)ei + (r3iVi - ri7V3)e2 
+ {riN2 - r2Ni)e3]. (AlO) 



where K{0, e) is a function of e and O (see the lower panel 
in Fig. lXTIl . 

Thus, unpolarized radiation produces only the third 
component of RATs, i.e. Q^s for the spheroid, while two 
first components Qei,Q<i2 vanish. 



APPENDIX B: MIRROR ON A POLE MODEL 
Bl RATs calculations 

Consider now a grain consisting of a square refiective mirror 
of side I2, attached to the spheroid by a pole of the length li 
(see Fig. [2]). Here we calculate RATs acting on the mirror. 
Its orientation in the grain coordinate system, ai,a2,a3, is 
characterized by a normal unit vector N, given by 



N = riiai + 7128.2, 



(Bl) 
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where ni = sin a, n2 = cos a with a is the angle between 
N and £12. (see Fig. [2]). The use of the grain coordinate 
system is appropriate as the spheroidal body determines the 
grain inertia. 

Since li ^ I2, the radius vector determining the position 
of each reflecting event on the mirror surface r is nearly 
parallel to aa, thus 



r = has. 



(B2) 



The orientation of the grain in the lab coordinate system is 
determined by O and /3 as follows, 

ai — cos Oei + sin Be2, (B3) 
a.2 ~ cos sin 0ei + cos Oe2] + sin jSes, (B4) 
aa = — sin /3[— sin Oei + cos Be2] + cos jSeg. (B5) 



Plugging equations (|B3IB5P into equations (|B1|I and (|B2|) 

we get 

N — ei(ni cos B — 712 cos (3 sin O) 
+ 62 (ni sin O + W2 cos /? cos Q) 

+ ean2 sin /3, (B6) 
r — li{ei sin (3 sin O — §2 sin /3 cos B + ea cos f3). (B7) 

Hence, RAT produced by the photon beam acting on the 
mirror is 

dVrad = 7r X AP = (-2pph)7FdA(k.N)[r x N]. (B8) 

Integrating over the full mirror, RAT becomes 

Trad = -2p'yFAi_ (k.N) [r X N] , (B9) 

where A± = j4|ei.N| is the cross section of the mirror with 
respect to the photon flux. 

Substituting pph, F, k = ei, and A± — i||ei.N| into 
equation (|B9[) . we get 



Trad = ^"''""^^^ (-|)iei.N|(ei.N)[(r2jVa - raiV2)ei 

+ (raiVi - riiV3)e2 + (riiV2 - r2Ai'i)e3]. (BIO) 

From equations (|B6P and (|B7[) . we obtain 



r2A^3 - 


rsNi 


h 






riiVa 


h 




riNi - 


r2Ni 



= — n2 sin f3 cos B sin f3 

— cos /3(ni sin B + 712 cos /3 cos B) 

= (— ni cos /9 sin B — n2 cos B), (Bll) 

= cos /3(ni cos B — n2 cos /3 sin B) 

— sin fi sin Bn2 sin /3 

= m cos /3 cos 6 — n2 sin 6, (B12) 

= m sin /3. (B13) 



Therefore, we can write 



T rad — 



■yUrad^^ll 



(Qeiei+Qe2e2 + Qeae3), (B14) 



where the RAT components are given by 
4/1 

Qei = r- |ei .N| (ni cos f3 sin B + 712 cos 6) [— ni cos B 

A 

+ (n2 cos /?) sin B], (B15) 
4/1 . - 

Qe2 ~ — |ei.N|(ni cos /3 cos 6 — n2 sin 6)[ni cos 6 

A 

— n2 cos (3 sin B], (B16) 



4/1 . - 

]e3 = r-|ei.N|ni sin /3[ni cos B — n2 cos f3 sin B]. 

A 



(BIT) 



Simplifying the above equations, we get 

Qei = -^1^1 cos B — 712 sin B cos f3\[nin2 cos^B 
A 



+ ^ cos /3 sin 2B - cos /3 sin 2B 



' ni722 sin^B cos^ 



(B18) 



4^1 2 2 

Qe2 ~ — 7-^1 cos B — 712 sin B COS /3| [jii cos /3 cos B 



nin2 2r, ■ r,r\ ^^"^'^ 
— - — COS p sm 2B 



— sm 2B + 712 COS p sm BJ , 
(B19) 

^ea = cos B — 712 sin B cos f3\ni sin /3[ni cos B 

A 



— 712 cos (3 sin 6] 



(B20) 



Averaging over the rotational angle (3 in the range 
[0, 2-k], we get 

Qei = ^^^iJi^(3cos^B - 1)/(B, a), (B21) 



47rZinin2 . „^ ^ 
Qe2 = r sm 2Bg(B,Q) 



(B22) 
(B23) 



Qe3 = 0, 

where f{&,a),g{Q,a) are fltting functions depending on a 
and B. The dependence on B characterizes the influence of 
variation of the mirror cross section on RATs. We will flnd 
these fltting functions in the following section. 



B2 Fitting functions 

As we have seen above, RATs can be decomposed into an- 
alytical terms and fitting functions which are functions of 
both 6 and a (see equations IB18I and IB19|) . Below we dis- 
cuss an analytical approximation to the fitting functions 
/(B,a),p(6,a). 

In the vicinity of B = 0,7r, |7ii cos B— 712 sin B cos /3| ~ 
rii, so this factor does not make /3 -averaging more involved, 
however. 

As 6 ~ 7r/2, equations (|B15P and (|B16|) can be written 

as 

Qei = — -!-| — 712 cos [3\[nin2 cos^B + — cos f3 sin 2B 
A 2 
2 

^ cos (3 sin 2B — 771772 sin^B cos^/3], (B24) 

^ 4^1 I air 2 2r\ ^1772 2 r, ■ r,, 

Qf,2 ~ T- \ — 7i2 cos pj[77i cos [3 COS B COS (3 sm 2( 

A 2 



711772 . „^ , 2 o . 2r-,l 

— - — sm 2B + 772 cos p sm BJ, 



(B25) 



Integrating equations (|B24|) and (|B25|) over /3 in a range 
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Figure Bl. Tabulated and fitting functions / for a = 45". 
Sliaded lines show the vicinity of the singular line cos = 1/ 



a = 45 







1 .0 - 
^, 0.8 - 

51 - 

0.6 - 

OA - 

0.2 7 

0.0 
-1.0 



tabulated function g 
fitting function g 



-0.5 



0.0 

COS0 



Figure B2. Tabulated and fitting functions g for a = 45". 
[0, 2-k], we get 

Ah 4nin2|n2| , 2^ o^ mofi\ 
Qei = — TT — ^ ~ 2), (B26) 



A 3 

ill 10nin2\n2\ 



sin 29. 



(B27) 



A 3 

Comparing equations (|B18|l and (|B19|l with (|B26|) and 
(|B27|) . we have 



/^/2(6,a) = 
5,^/2(6,0) = 



|cosa|(5 cos^G — 2) 
37r(3 cos^e - 1) '' 
10|cosq| 
3^ ' 



(B28) 
(B29) 



The fitting functions /,r/2(0, q) and the tabulated func- 
tion f(<d,a) for a — 45° are shown in Fig. ([Bill . 

Interestingly enough, the approximation of the fitting 
functions obtained in the vicinity of ~ 7r/2 presents a 
reasonable approximation for the fitting function / over the 
entire range of O. Indeed, Fig. IBll shows the tabulated func- 
tion a) obtained by substituting Qei, Qe2 resulted from 
numerically averaging equations HB15|I and (|B16|I into equa- 
tions (|B18|) and (|B19|) fsolid lines) and the fitting function 
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Figure B3. Fitting function / for two tilted angles of the mir- 
ror in the body system: upper and lower panel correspond to 
a = 30" , 60" . Shaded lines show the vicinity of the singular line 
I cos e| = 1/^3 



/7r/2(0)Q^) given by equation HB28|I (dashed line) for a par- 
ticular case of a = 45" . It is shown that the fitting function 
/ir/2(05«) has a good agreement with f{Q,a). Similar re- 
sults are also found for other a angles (see Fig. IB3|) . 

However, the fitting function given by equa- 

tion (IB29|) is independent of Q and therefore does not pro- 
vide a good fit with the tabulated function g. In this case, 
we perform a quadratic fitting for tabulated functions for 



30", 45" and 6O" are 



different a angles. The results for a - 
given by 

g{e, a = 30°) = 1.4675 - 0.5 cos^e, (B30) 
g{e, a = 45°) = 1.191 4- 0.1382 cos^S, (B31) 
g{e, a = 60°) = 0.864 + 0.869 cos^O. (B32) 

Tabulated and fitting functions g are shown in Fig |B2l for 
a = 45° (see Fig. IB4l for g corresponding to other angles a). 

Figs (|B3l) and (lB4| show f,g for a = 30" and 60°. Fig. 
IB3I shows that there is a change in sign of the tabulated 
function / as | cos e| l/VS between a = 30° and 6O". 
This stems from the fact that, for a < 45", the influence of 
A± shifts the zero of Qei from | cos Q\ — l/\/3 to | cos 0| < 
l/VS. As a result, when the term 3cos^0 — 1 approaches the 
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Figure B4. Fitting and tabulated functions g for two tilted an- 
gles of the mirror in the body system; upper and lower panel 
correspond to a = 30", 60". 



zero from the left of cos Q = — 1/\/3, for instance, then Qei 
still positive, which causes / rises sharply to the positive 
direction as in Fig. IB3b pper. 

As an example, in the paper we mostly discuss the case 
in which the mirror is tilted by an angle of a = 7r/4 for 
which there is the best fit between the fitting function /.^/2 
and the tabulated function / except a narrow range in the 
vicinity of | cos e| = 1/73 (see Fig. IBip . Similarly, g given 
by equation (|B31|I also provides a good fit for Qe2 (see Fig. 
B2[l . Therefore, RATs acting upon the mirror can be roughly 
approximated by equations (|B18|I - (|B19|I with f,g are given 
in equations HB28|I and (1B31|) . 

Needless to say, that such an approximation entails er- 
rors in calculating the torques Qei and Qe2- Fig. IB5I shows 
the torques with /, g being the tabulated and fitting func- 
tions. It can be seen that there is a very good fit for the 
component Qe2, but there exhibit some deviation for the 
component Qei toward the range cos G ~ ±1. These errors 
could be decreased by using more sophisticated analytical 
fits, e.g. piecewise analytical fit. However, this would de- 
crease the heuristic value of the formulae. Therefore while 
using the approximate fits (|B28P and (|B31|) for some ana- 
lytical calculations below, for most of the quantitative esti- 
mates, including those related to evolving the phase trajec- 
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Figure B5. Fig shows the comparison of RAT components. 
Dashed lines show RATs obtained from the analytical approxima- 
tion given by equations | |B18| |- ||B19| I and solid lines show RATs 
obtained by numerically averaging equations IIB15t - IIB16l l over /3. 



tories, we tabulate the fitting functions / and g. Naturally, 
the latter is equivalent to the direct use of averaged equa- 
tions (|B15|) - pi6|) for RAT components. 

What is the purpose of having the analytical form of 
the fitting functions, if they do not provide exact fits for 
the torque components (see Fig. IB5|l ? We may see the ac- 
tual RATs acting on irregular grains (e.g. Fig. [5]) also differ 
somewhat from the corresponding AMO RATs. These dif- 
ferences, however, do not change substantially, as we show 
in the rest of the paper, the alignment of irregular grains in 
comparison with the alignment predicted by AMO. Thus we 
believe that the slightly distorted AMO with the appropri- 
ate analytically fitting functions should also reflect correctly 
the generic properties of the RAT alignment. 



B3 Dependences on a 



Fig. IB6I shows that the variation of the mirror orientation 
can change the ratio of their maximum QS^^/Q™^- Here 
we denote Q™^ the maximum of IQei] which is exactly its 
amplitude at O = or tt, and QS"^ the maximum of Qe2- 
In Fig. |B6]it is obviously shown that Q^'^/QS"'" increases 
with a increasing due to the increase of the cross-section A± . 
However, the ratio Q]PC^ /Qe2 is limited to the upper limit 
of 1.3. We will see later that the ratio plays an important 
role on the alignment that is discussed in § 17.51 

Moreover, the variation of a does not only give rise to 
the change of the ratio QTi"^ /Qlh"'^ , but it can also affect 
the symmetric functional form of RATs. Our calculations 
show that for a smaller than 20" or larger than SO", the 
RATs have some changes in the functional form of the Qei 
and Qe2 compared to what is shown in Fig. [3] 

Let us explain why for a small, the symmetry of RATs 
for AMO is affected. Consider equation (|B18|) with ni — 
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where rj is radius of dipole, Ej , Fj are electric field at the 
location of, and radiative force which acts on the j*'' dipole. 
Radiative force F present in equation (|C1|) is produced by 
the gradient of electric field in the grain and the Lorentz 
force due to vibration of electric dipole in magnetic field. 

Radiation field inside the grain consists of that of inci- 
dent and scattered light. Each dipole receives the incident 
light which induces its vibration and scattered light pro- 
duced by all electric dipoles except the dipole under study. 
Hence, the total RAT can be written in a different form 



ratio determines the existence of high attractor points and the 
shift of a crossover point as we describe in i]6.6l Shaded areas show 
the range of a in which the symmetry of RATs are significantly 
affected by the cross section. 



Frad ^ inc ~t- F^ca- 

RAT efficiency Qr is defined as foUowings, 



Trad — 



Q 



sm a, 712 = cos a 
'ill 



A 



+ 



cos f3 sin 20 



2 ' 2 

sin a cos a sin^O cos^/9]. 



cos P sin 20 



(B33) 



The first term is the cross-section A±, and the second term is 
Qei without cross-section. The second term when averaged 
gives rise to Qei ~ 5 cos^B — 2. It indicates that if the cross- 
section is constant, Qei is fully symmetric. Consider a small, 
cosa is significant, and equation (|B33|I reduces to 



cos Q sin O cos /3j[ sin a cos a cos O 



H — cos f3 sin 20 

— sin Q cos a sin^B cos^/3] 



cos /3 sin 2B 



(B34) 



Obviously, when averaging over l3, the absolute term con- 
taining sin B contributes to modify substantially the sym- 
metry of the resulting torques. The same problem also oc- 
curs for Qe2- We found that for both a ^ 20" or q ^ 80", 
the functional forms of RATs are very influenced. Ifowever, 
for a within [20^,80"], their functional forms are not much 
different with what shown in Fig. |3] for a — 45". 



APPENDIX C: RATS WITH DDSCAT 

In the discrete electric dipole approximation (Draine & 
Flateau 1994), a grain is presented as an ensemble of electric 
dipoles. The interaction between the electric field of incident 
light and the dipoles produces radiative forces and torques. 
RATs produced by radiation on a grain consisting of A'' elec- 
tric dipoles are 



(CI) 



-7Qr, 



(C2) 



(C3) 



where 7 is the anisotropy, and ux is the energy density of 
radiation field of wavelength A. Hence, equation (|C2|I can be 
rewritten 



Qr — Qinc -t- Qsca. 

Here Qi„c = Qabs and Q^ca are given by (DW96) 

AT 

4fc 



(C4) 



sin a cos B — cos a sin B cos f3\[ sin a cos a cos B 



Re^Pj(O) X Ei„c,oe" 

N 

-ik X y^rj[pj.Eine,o]e''' 



Qsca — 



2— / d^ReiSEVB +S'bVe) 

ef f i^incO I J 



(C5) 



(C6) 



In above equations, Ei„c denotes electric field of incident 
light, Se,Sb, Ve,Vb are given by 



2i 

-^n]-Pj(0)ea;p(-ifcn.rj), 



Sb ~ r. ^^Pj(0) X rjexp{—ikn.rj 



Ve 



Vb = 



^Pj(O) - n[n.pj{0)exp{-ikn.rj) 
-N X Ve, 



(C7) 

(C8) 

(C9) 
(CIO) 



where pj is electric dipole moment, k is wave number, r 
is the radius vector and n is normal unit vector. 

We use DDSCAT code to compute separately the com- 
ponents arising from absorption, scattering and total RAT 
for a grain in which its direction with respect to k is de- 
termined by the angles B,/3, $ (see Fig. Here B is the 
angle between ai and k; /3 is the rotational angle of the grain 
around ai; and "I> is the precession angle of ai about k. For 
our study (see Table [TJ, we compute RATs for the spec- 
trum of the interstellar radiation field, over 21 directions of 
6 from to TT and 21 values of (3 from to 27r, at $ = 0. 

RAT can be decomposed into components in the scat- 
tering system via 



Qr = QelBl + Qe2e2 + QeiSi, 

where ei, 62, 63 are shown in Fig. [1] 



(Cll) 



i=i 



Mean radiative torque over wavelengths, Q(6, /3, $) is 
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defined by 

— _ / QxUxdX 
J u\dX 



(C12) 



Since (3 varies very fast due to the swift rotation of the 
grain around the axis of major inertia ai, we can average 
RATs over /3 from to 27r. 



APPENDIX D: ATTRACTOR AND REPELLOR 
POINTS 

Here we derive the condition for which a stationary point 
becomes an attractor and a repellor point foUowing the ap- 
proach in DW97. Assuming that the stationary point has the 
angle with respect to magnetic field and angular momen- 
tum Js, one can expand the right hand sides of equations 
of motion around this point. As a result, equations (|66|l and 
((67)) give 



dt 



Js 



(J-Js), (Dl) 



{^-^s)-{J-Js), (D2) 



where {H{(,)}4,, {F"{£,)}<i> ^-re spinning and aligning torques al- 
ready averaged over the precession angle (j). Since for station- 
ary points ^s,Js, we have {F(5s))</. = 0, {H{^s)),p — Js = 0, 
equations (|D1[) and (|D2I) become 



dt 
dui 



dt 
where 



= A{i-is) + B{j-js), 

= C{^-^s)+D{J-Js), 



B 



JdC 

JI ' 

_ d{H)^ 
C - -^{is,Js), 

D = -1. 

To have an attractor point, one requires 

A + D <Q, 
BC~ AD < 0. 

In other words, 

d{F)^ 



{H)^d^ 



{^s,Js) < 1, 



H 



d{F), 



{is,Js) < {F)^^^{^s,Js) 



d^ " ^ ' dC 

where we have substituted Js ~ H{^s)- 



(D3) 
(D4) 

(D5) 

(D6) 

(D7) 
(D8) 

(D9) 
(DIO) 

(Dll) 
(D12) 



APPENDIX E: EFFECTIVE GRAIN SIZE FOR 
AMO 

For phase maps, we numerically average equations (|10|l - p2| ) 
to obtain exact RATs for AMO, rather than using approxi- 
mate formulae as in the analysis. However, for this case, the 



absolute magnitude of torques matters. Therefore, we nor- 
malize AMO in the following way. Assuming that the size I2 
of the mirror and h are chosen so that the RAT for AMO 
has the magnitude equal to that of irregular grain of an ef- 
fective size aeff. Thus, following equations Q and ((9)1 we 
have 



Qamo 



XUradoiff 



Qddscat, 



(El) 



2 — - 2 

where Qamo and Qddscat are the magnitudes of RATs for 
AMO and an irregular grain, respectively. We can simplify 
further by normalizing RAT components over the maximum 
of Qei, and let Qamo, Qddscat equal the maximum of 
Qei for AMO and irregular grain. As a result, for AMO, we 
have Qamo ~ IQei]™"^ = -'-^^"^"a'l ^ g^j^^j substituting into 
equation HE1[I . we get 

^ ' " (E2) 



'DDSCAT- 



We can define the effective size of AMO subject to the mirror 
size and the rod length as 



2 



Si'i). 



(E3) 



Equations (|E2p and HE3|) enables us to find the effective size 
of AMO that produce the same the RAT magnitude with 
the irregular grain of size fle//, given by 



2 

O-AMO 



1 



deff-QoDSCAT- 



(E4) 



As we discuss in ijlOl that the magnitude of RATs for 
irregular grains can be crudely approximated as 

A 



= 0.4( — )" 

a-eff 



(E5) 

3 or 4 for A > 1.8 cte// 



where r] = for A < 1.8 fle// and ri - 
(see 

Equations HE5[) and (IE4p allow us to roughly estimate 
the the effective size of AMO, aAMO as a function of the ef- 
fective size of irregular grain, a^f / , provided that the wave- 
length of radiation field is known. 
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